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Abstract 



In this master thesis I discuss how the weak three-body decay tx^-k^K^ is 

related to the strong interactions between the final-state particles. Current issues in hadron 
spectroscopy define a region of interest in a Dalitz plot of this decay. This region is roughly 
characterized by energies from to 1.6 GeV of the -k^-k^ and the ti'K^ system, and 
energies from 3 to 5 GeV of the n^K^ system. Therefore I propose to use for the former 
two systems elastic unitarity (Watson's theorem) as an approximate ansatz, and a non- 
vanishing forward-peaked amplitude for the latter system. I introduce most basic concepts 
in detail with the intention that this be of use to non-experts interested in this topic. 
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Chapter 1 



Introduction 



1.1 Dalitz plots and hadron spectroscopy 

Dalitz plots are representations of a three-body decay, 

X abc, (1.1) 

in a two-dimensional plot. The two axis of the plot are nowadays usually the invariant 
masses squared of two of the three possible particle pairs, i. e. for instance^ 



Sab = {Pa + Pb) 
Sac = {Pa + Pc) 



(1.2) 



One can also choose for the axis the invariant masses not squared or, as it was done mostly 
in the first Dalitz plots, the kinetic energies of two of the three decay products, for instance 
Tfe and T^. These coordinates are equivalent in the sense that they are linearly related, e. g. 

Sab = rn\+ml-2mx{mc + T^). (1.3) 

Dalitz plots owe their name to Richard Dalitz who developed this representation tech- 
nique in order to analyze the decay Tr+Tr+Tr" jUEj.^ The basic idea is to assign each 
decay event its coordinates with respect to the two axis of the plot. One thus can obtain 
a "landscape" where the "mountains" correspond to a lot of events and the "valleys" to 
very few or no events. 

1.1.1 "Old" and "new" Dalitz plots 

Originally, the main application of Dalitz plots was to use it in order to determine spin 
and parity of the decaying particle. A prominent example is the decay 

TT+vrV-, (1.4) 



^Somc of the kaons were then caUed "r-meson". See also P- 141]. 
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see figure ll.ll Spin and parity of tlie decaying particle can be read off from cliaracteristic 
patterns in tlie Dalitz plot as discussed prominently in ref . jl] , see figure 11.21 

In recent years Dalitz plots came to the fore again mainly in studies of D and B decays, 
see for instance [7j and jS]. These "new" Dalitz plots are distinguished from older Dalitz 
plots like the one in figure 12.11 by that they contain a number of events of 1000 to 2000, 
which reduces statistical fluctuations. Such modern Dalitz plots (figure II. 3p are therefore 
sensitive to the presence or absence of slight variations in the event distribution. Older 
Dalitz plots like the one in figure 12.11 could be used rather for "resonance hunting" . One 
searched for bands in the plot which can be associated to resonances. 

Thanks to the sensitivity of the new Dalitz plots, they provide new input for problems 
in hadronic spectroscopy, see for instance P and [TU] . 

1.1.2 Issues in the spectroscopy of hadrons 

Today's controversial issues in hadron spectroscopy include the isoscalar scalars, more 
precisely the states with the quantum numbers of the vacuum, 

J^^ = 0++. (1.5) 

One reason why these states are of particular interest is that two-pion exchange, which 
has J'^*" = O^"*", is the next-to-longest range contribution for the strong force p. 1185]. 
The longest range contribution is the exchange of pions, the lightest hadrons. The isoscalar 
scalar sector is also relevant for gluehalls. Various models (see references in ^2j) predict 
the lightest glueball to have the quantum numbers of the vacuum. 

To learn more about the scalars one investigates the component with total angular 
momentum J = (5" wave) and isospin 7 = (iso-scalar) mainly of the vrvr scattering 
amplitude, see figure I3.fj| top panel. For the low energy region from zero to about 500 
MeV there is Chiral Perturbation Theory that can state various results. However, the 
chiral expansion is hardly valid for energies up to about 1.6 GeV, which are to be taken 
into account in discussing the scalars. 

The characteristic feature of vrvr isoscalar S wave cross section, called red dragon by 
Minkowski and Ochs, can be explained by the negative interference of the /o(980) and the 
/o(1500) with the broad glueball with J^*-^ = 0"'"+, ^21- Negative interference means that 
the rapid phase shift of it does not start from zero as for a pure Breit-Wigner resonance but 
from a background phase of about already it/ 2 and then again from 37r/2, see figure ITTI 
on page |321 Alternative explanations postulate the a (also known as /o(600)) and the 
/o(1370). The two alternatives yield at least to different scalar nonets as described in pj: 

• A scalar nonet with /o(980) and no lighter particles, including ao(980), J<'o(1430) and 
/o(1500), and 

• a scalar nonet with /o(980) and no heavier particles, including ao(980), k(850) and 
(t(600). 
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Figure 1.1: Dalitz plot for the decay ^ tt'^it^'it , from [S] cited in p. 320]. 

1.1.3 Weakly coupled decays and strong interactions 

Remarkably enough, vrvr scattering cannot be studied experimentally by scattering a pion 
beam off a pion target. This is impossible because of the short lifetime of about 10~^s 
for the charged pions and 10~^^s for the neutral pion. A pion target would therefore have 
disappeared even before any scattering can happen. Chew and Low ^2] proposed around 
1960 a surrogate laboratory for unstable targets. Using basically the idea of Chew and 
Low a big part of the experimental input for irn scattering comes from experiments such 
as np TTirn, see also section IH. 2. 31 

Another surrogate laboratory are three body decays. A prominent example is the decay 
J/V' — 4>7T7T. This is a process of the strong interaction but it is weakly coupled.^ More 
recently three-body decays by the weak interaction of B and D mesons have attracted 
attention. Dalitz plots from these decays that were produced in the first place for studies 
of CP violation serve with their high statistics as another source of a priori very precise 
experimental input. To use such Dalitz plots for the study of hadronic two-body scattering 
one has to relate the weak amplitude of the three-body decay to these hadronic amplitudes. 
To profit from the high statistics and to extract information not only from clearly visible 
peaks and dips it is important that the resulting ansatz for fitting the Dalitz plots is 
sensitive also to slight variations in the event distribution. 

The ansatz that I will try to elaborate is guided by unitarity in two respects. For 
the amplitudes where we are interested in the energy region where at least approximately 

■^Thc Particle Data Group reports a branching fraction of 8 x 10^'^ for the mode ^tt+tt" |14j . 
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Fig. 2. Regions of the Dalitz plot where the density must 
vanish because of symmetry requirements are shown in black. The 
vanishing is of higher order (stronger) where black lines and dots 
overlap. In each isospin and parity state, the pattern for a spin of 
even integer is identical to the pattern for spin /, provided 
/^2, (Exception: vanishing at the center is not required for 



Figure 1.2: Characteristic patterns in a Dalitz plot from which spin and parity of the 
decaying particle can be determined; from |3]. 
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only the elastic channel is open, unitarity leads to a form of the decay amplitude that 
has the same phase as the strong interaction amplitude of the hadrons in the final state. 
This result is known as Watson's theorem ^3]. For the amplitude in higher energy regions 
inelastic channels have to be considered and the characteristic high energy behavior of the 
elastic amplitude may explain the main features. 1 will discuss whether for the high energy 
amplitude it is again unitarity that can be used to obtain via the optical theorem a rough 
idea of its behavior. 

1.2 Outline 

In chapter 2 I will introduce basic features of Dalitz plots. I will briefly discuss some 
points concerning background suppression. Then I will stress the possibility of interference 
between amplitudes and how in the presence of interference branching fractions are to be 
understood. 

In chapter 3 the S and T matrix are defined. I discuss normalization and completeness 
of particle states, give formulations of unitarity conditions and derive the optical theorem. 
The detailed account of normalization allows me also to derive an equation defining the 
boundary of a Dalitz plot. 

Chapter 4 deals with the kinematics of a three-body decay. There again one obtains 
equations for the boundary of a Dalitz plot. This time not through phase space consid- 
erations but through the definition of two-body subsystems in the three-body final state. 
These two-body subsystems are of interest because they can be related to the strong 
two-body scattering amplitudes in the final state and can be used to define alternative 
coordinate systems in a Dalitz plot appropriate for a partial-wave analysis. 

In chapter 5 I will elaborate on the relation between the decay amplitude and the 
scattering amplitudes of the final state particles. For the scattering amplitudes I will first 
apply the constraints from unitarity (Watson's theorem and the optical theorem). Then I 
try to see what consequences the known features of the (strong) two-body amplitudes have 
on an ansatz for the (weak) amplitude of the three-body decay. 



Chapter 2 



Background, resonances and 
interference 

2.1 Background suppression 

Consider as an example the study of the charmless K^n^ii^ decay JB]- First, one 

has to detect the three particles with an appropriate detector. But not all i^'s and vr's that 
are detected during the experiment may be produced in a charmless B decay, and some 
particles may be misidentified as i^'s and vr's. 

The K^ir^n' final states that come from other reactions than a B decay can be vetoed 
by requiring that the total invariant mass of the three particles differ not much from the 
invariant mass of a B^. By a similar criterion, then, one can further exclude the final states 
which are produced through a charmed (instead of charmless) B decay (e. g. B^ £)°7r+ 
followed by D° ^ A'+vr") p. 7]. 

2.1.1 Sidebands 

The final states where the invariant masses are in a certain neighborhood of the B invari- 
ant mass belong to the signal region; other K~^7r~^7T~ final states belong to the so-called 
sidebands. The events in the sidebands are background events. But not all events in the 
signal region are signal events, i. e. they do not all come from a charmless B~^ decay. In 
other words, not all background events are excluded by a sideband criterion. In particular, 
the background that comes from misidentification is still present. This sort of background 
is usually estimated by an extrapolation from the sidebands: When one has excluded all 
events from the sidebands region, the question remains which events in the signal region 
do still not come from a charmless B~^ K^ti^tt^ decay. To answer this question one 
may try the following procedure: Assume that the background that is not excluded by 
a restriction to the signal region comes only from misidentification. Try to estimate how 
much misidentified events there are in the sidebands. Assume that the same number of 
misidentified events is also present in the signal region. Subtract this number from the 
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total events in the signal region and you will be left with events that indeed come from a 
charmless decay K^ti^-k^ . 

2.1.2 Background and correlations 

A Dalitz plot analysis consists essentially in extracting information from deviations from 
homogeneous event distribution over the Dalitz plot. A band of higher than average density 
of events parallel to one of the axes, for instance, is the correlation between the four- 
momenta of two out of the three final state particles. The information one can extract 
under certain circumstances is that the two particles in question form a resonance. Another 
type of information extraction from correlations in this context is partial wave analysis. 

In order to extract information from correlations one has to subtract the background in 
a way that does not give rise to correlations that are not due to the physics of the process. 
To avoid such spurious correlations one has to take care of two points in particular: 

Partial wave decomposition of the background. Also the background amplitudes 
should be decomposed into partial waves. This is often not done. As a positive example 
in this respect I can cite ref. ^Zj, where the function that characterizes the acceptance of 
the detector is expanded in spherical harmonics. 

Forming triplets of signal events. A point in a Dalitz plot represents a triplet of 
signal pions and kaons that come from the same B. Assume that constraints like energy 
and momentum conservation or detection time do not determine which pions and kaons 
make up a Dalitz plot event. Assume as a simple example that there is a set of 8 7r~'s, 
a set of 7 -7r"'''s and a set of 9 K^''s that are candidates for coming from the three-body 
decay of one of 5 -B^'s, see table ITTl This means that the background estimation gives 
a ratio of signal to background of 5/8, 5/7 and 5/9, respectively. Due to the possible 
quantum mechanical superposition of background and signal amplitudes it is impossible 
to say of an individual event whether it is a background or signal event and — if a signal 
event — from which of the 5 -B's it has come. But in order to avoid spurious correlations it 
is not sufficient just to fix the signal to background ratio of candidates for forming a Dalitz 
plot triplet. One must not form, for instance, pairs of pions for the ti'tc^ mass projection 
and 7r~K'^ pairs for the 7t~K~^ mass projection where the and the vr"*" are not assigned 
to belong to the same Dalitz plot triplet. As an illustration see table fT/H 

2.1.3 Experimental difficulties 

Particular experimental difficulties are encountered in obtaining a Dalitz plot in decays 
where one of the three decay products is a photon or where some decay products can 
further decay into photons, as in 



J/ijj 'jTTir and B — > Krjir. 



(2.1) 
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Table 2.1: An example of forming pairs for mass projections and Dalitz plot triplets that 
is consistent with a fixed signal to background ratio but can lead to spurious correlations. 



The J/ip can also decay via a p and a vr^ into a four-particle final state tttt^^, 

J/?/; -> p7r° ^ (7r7r)(77). (2.2) 

Photons escape quite easily any detection. If this happens the reaction J/ip ^ pn^ 
(7r7r)(77) is erroneously identified as a Dalitz plot event J/ip ^ 'jim; cf. jTTJ p. 1193], 
where these experimental difficulties are briefly recalled, and the work of the Mark III 
collaboration (e.g. ref. [HI)) one of the most successful in overcoming these difficulties. 

In the case oi B ^ K-qir, on the other hand, the loss of a 7 may have as a consequence 
that a Dalitz plot event is not identified as such. This is because one way to detect the 
final state 77 is via its subsequent decay into two photons. 



2.2 Dalitz plots 
2.2.1 Basic features 

Dalitz plots have been introduced as a convenient method of analyzing reactions with three- 
body final-states in ref. jl], see also ref. Consider a three-particle final-state consisting 
of the three particles abc. Out of these three particles one has three possibilities to form a 
pair of two particles: ab, ac and be. For a Dalitz plot the invariant masses squared of two of 
the three pairs are represented as the two coordinate axes.^ Each of the three-body final- 
states produced in an experiment can be represented as a dot in the plane defined by these 
coordinate axes. Theoretically such dots are confined to a certain boundary depending on 
the kinematics of the reaction because outside this boundary the phase space volume is 
zero. For the interior of this boundary the phase space volume is constant, see section 
Therefore, the density of the points inside the boundary is a constant multiple of the 

^Originally, see e.g. ref. each coordinate axis represented the kinetic energies of one of the three 
particles. 
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reaction matrix element squared. Thus, every structure in the density of the plots is due 
to dynamical characteristics of the reactions and not of kinematical origin. 
If the reaction proceeds via an almost stable two-body intermediate state, 

X {ab)c abc, (2.3) 

the three-body final states will be such that the distribution of the invariant mass of the 
pair ab is centered around the invariant mass of the two-body intermediate state. This 
manifests itself in the Dalitz plot as a band of higher than average density of points, see 
figure 12.11 

A three-body decay X — > abc has three two-body channels: ab, ac, and be. If e. g. the 
pair ab can form a quasi-stable state, a resonance, this appears in the Dalitz plot as a band 
along the region of the plot where the the two-particle invariant mass of the pair ab is equal 
to the mass of the resonance.^ If e.g. the pair ab may form a resonance, and one of the 
axes of the Dalitz plot is the invariant mass of the pair ab, then the band is perpendicular 
to this axis. 

The same type of particle, say ab, may form more than one resonances. The various 
resonances then appear in the simplest case as parallel bands of a certain width in the Dalitz 
plot at respective positions given by the mass of the resonance. In general, interference 
effects can generate different structures. 



2.2.2 Interference of a resonance with itself 

Consider as an example of an interference effect a three-body decay where two of the three 
final-state particles are of the same type, 

X abb. (2.4) 

Suppose further that a particle of type a and a particle of type b may form an almost stable 
intermediate state, a resonance R*, say. Moreover, the channel via the resonance R* be 
the only process by which the particle of type X decays into a particle of type a and in two 
particles of type b. Then the decay amplitude is a superposition of two amplitudes: The 
intermediate state can be made up of the particle of type a and one or the other particle 
of type b. 



X 



f 


a 




b 


b 1 




a 


b 


b 


1 



abb. 



(2.5) 



In the Dalitz plot of such a three-body decay, with mat of one combination ab on one 
axis and m^f, of the other combination on the other axis, there are two resonance bands 



•^This statement needs qualification, see e.g. ref. |2()| . 
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Figure 2.1: Bands in a Dalitz plot oi K p ^ tt^Att^ due to the pion-lambda resonance 
S±(1385); from Also shown: Histo grams of projections on the axis representing 

kinetic energies of two out of the three final-state particles. 
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perpendicular to one of the axes, respectively, at the position given by the mass of the 
resonance R*. Depending on the kinematics of the process the two bands may or may not 
overlap. In case they do, there is constructive interference of the resonance amplitude in 
one of the two resonant two-body channels with the amplitude in the other. Thus in the 
region of the overlap the total intensity is not twice the intensity of an isolated band but 
four times because it is the amplitudes and not the intensities (i. e. the amplitudes squared) 
that are added. Two examples with even three overlapping bands are shown in figures 12.21 
and 12.31 In figure 12.21 the bands are regions of higher than average density, i. e. peaks. In 
figure the light blue bands due to the resonance /o(980) are regions with lower than 
average density, i. e. dips or valleys. Besides the overlap of bands this figure thus shows 
remarkably that resonances (here the /o(980)) do not necessarily show up as peaks in a 
cross section. 

2.3 Branching fractions 

For a lot of unstable particles there is more than one possible decay product. A particle 
X may e. g. decay in two particles ab, or cd, or in three particles efg. The probability for 
the particle X to decay in one of the possible products is the respective branching fraction. 
The most likely decay could for example he X ^ ab with a probability of 80%. With 15% 
probability the particle X may decay in the particle pair cd, and with 5% probability into 
efg. If there is only one possibility for the particle to decay, the branching fraction for this 
decay is 100% and for all other conceivable decays 0% — or one may prefer not to speak of 
branching fractions at all in such a case. 

Dalitz plots can be used to determine branching fractions for three-body decays with 
almost stable two-body intermediate states. Consider as an example the following three- 
body decay that can proceed via three different resonances: 



One may then ask what the respective probabilities are for these three possibilities. These 
probabilities are the (exclusive) branching fractions for the three processes respectively. 
The probability for the three-body decay X — > abc to occur at all, i. e. via any intermediate 
state, is the inclusive branching fraction for the process X abc. 

The task of determining branching fractions from a Dalitz plot may be seriously com- 
plicated by interference effects: 

[. . . ] we determine the exclusive branching fractions neglecting the effects of 
interference. The uncertainty due to possible interference between different in- 
termediate states is included in the final result as a model-dependent error. [. . . ] 
We find that the model-dependent errors associated with the wide resonances 
introduce significant uncertainties into the branching fraction determination. 




(2.6) 



El 



2.3 Branching fractions 
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Figure 2.2: Overlapping bands showing the interference of the p resonance with itself 
(different charge states) in pp tt^tt^'k^; from |25 cited in [22 P- 36]. 
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m2(jc0ji;0) [GeV2] 

Figure 2.3: Overlapping of valleys due to the resonance /o(980) in pp Stt^; from [SB]- 
Red and blue regions correspond respectively to high and low event densities. 
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2.3.1 Branching fractions, superposition and interference 

The basic idea of branching fractions is to give a quantitative answer to the following type 
of questions. "If various processes are possible under given circumstances: How likely are 
the respective processes to occur?" Or: "If various processes could have occurred: How 
likely is it that one particular process has occurred?" 

The quantum mechanical principle of superposition, however, implies that such types 
of questions do not always have an answer — except for the reply that the question is 
meaningless. This is for example the case in the notorious double slit experiment. If both 
slits are open as a matter of principle one cannot tell for a particular particle whether it 
has arrived at the screen via the upper or the lower slit. But still it seems that one can 
say something more about the intermediate state of a double slit experiment than just: 
"The particles arrived somehow at the screen." In some sense the two possibilities (passing 
through the upper or lower slit) are on an equal footing. So one would like to say that it 
is equally probable that the particle has taken the upper slit or the lower slit, but at the 
same time, because the particle was in a state of superposition of the two possibilities, it 
has not taken one or the other route nor both. How can one solve the dilemma between 
superposition and the basic idea of a branching fraction? How can one define branching 
fractions in a sensible way even in the presence of superposition? 

A possible solution is via an amplitude analysis. One can identify the intensity pattern 
as the square of a total amplitude which in turn decomposes into a sum of two amplitudes 
in the double slit case. One of these amplitudes represents the process via the upper slit 
the other amplitude represents the process via the lower slit. The respective squares of 
these two amplitudes are equal. It is in this sense that the two processes are on an equal 
footing. In this sense both processes (passing through the upper or lower slit) are equally 
probable although it is not the case that one or the other process has occurred but rather 
a superposition of both. Because the amplitudes of both processes contribute equally 
to the total amplitude the two processes should be assigned equal branching fractions, 
without interpreting this assignment as "in one half of a certain number of repetitions 
of the experiment the particle passes through the upper slit and in the other half of the 
repetitions through the lower slit". Such an interpretation would be at odds with the 
superposition principle of (standard) quantum mechanics. 

Let ip be the total amplitude of the process of a particle reaching the screen in the 
double slit setup, and ipi and ?/'2 the amplitudes for passing through the upper or the lower 
slit, respectively, 

= ^1 + ^2- (2.7) 
IV'P is the total intensity, which we set as unity of the branching fractions, 

1 = = Itfji + = lifjil^ + I^AaP + 2Re(^AiV^*)- (2-8) 

If as in the double slit experiment the two amplitudes interfere, we have Re('?/;i'?/;2) 7^ 0, 
hence l^/'ip + |^/'2p 7^ 1- This shows that when there is interference, branching fractions do 
in general not add up to unity. 



16 



Chapter 2. Background, resonances and interference 



The basic idea of branching fractions is to assign probabihties to various possibihties 
under given circumstances. Due to quantum mechanical superposition the totahty of pos- 
sibihties is not a complete set of alternatives in the sense of ordinary language, where 
exactly one of the alternatives is realized. If, moreover, the superposed amplitudes inter- 
fere the sum of all branching fractions is not one, which shows again, that it is not complete 
alternatives in the ordinary sense that branching fractions are assigned to. 

2.3.2 Fit fractions in Dalitz plots 

In Dalitz plots interference between resonant amplitudes show up as overlapping bands. 
As discussed in section in the presence of interference the definition of a branching 

fraction is less obvious than it may seem at first sight. A very early determination of 
branching fractions from a Dalitz plot with two overlapping bands is ref. |25j . 

The method applied there is roughly the following. The mass projections of the Dalitz 
plot are fitted with a Breit-Wigner distribution without taking into consideration the events 
that fall in the overlap regions. The density of events in the resonance bands thus obtained 
is then extrapolated into the region of the overlap, assuming constant density of events 
in the resonance bands. The observed number of events in the overlap region is in good 
agreement with this extrapolation. This shows that the two overlapping resonances do 
not interfere. Therefore, the first fit can be corrected for the overlap regions assuming 
incoherent addition of the two amplitudes representing the overlapping bands. 

But what about cases with interference? As mentioned at the beginning of section 12. HI 
the effect of interference on branching fractions is not always negligible. I cite again from 
ref. 1211: 

We find that effects of interference between different two-body intermediate 
states can have significant influence on the observed two-particle mass spectra 
and a full amplitude analysis of the three-body B meson decays is required for 
a more complete understanding. This will be possible with increased statistics. 

It seems that recent experiments provide sufficient statistics to allow for such a full ampli- 
tude analysis [121 122] • 



Chapter 3 

Scattering amplitudes and particle 
states 

3.1 S- and T-matrix 



Let |Q!)in and |a)out be some appropriately defined states that are identical to the free- 
particle states \a) (without subscript) in the respective limit t —oo t +00. The 
amplitude for the initial state |a)in to be found as the final state |a)out is given by the 
scalar products between these states. All possible scalar products can be collected in to 
rows and columns of a matrix, usually called S (for scattering, I suppose). 



We can define an operator S such that sandwiched between free-particle states it yields 
the corresponding element of the S'-matrix: 



The amplitude for the transition of state to the state \a) is the superposition of the 
amplitude which represents the quantum mechanical collapse of the state in its 

component \a), and an amplitude that represents a transition caused by some interaction. 
This latter amplitude is described in terms of the elements of the T-matrix 



Sa/B = out (« 1/5); 



(3.1) 



Sap = {(y\S\P). 



(3.2) 



{a\T\P) 



(3.3) 



such that 



Sa/3 = {(y\S\P) 



{a\p) + i{a\T\p) 



(3.4) 



i. e. we define the T-matrix through the relation 



S =l + iT. 



(3.5) 
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3.2 Cross sections and decay rates 

In this section I try to establish — in some detail but without going into the very foundations— 
the relation between the T matrix elements and cross sections and decay rates. In partic- 
ular I will focus on a decay of type — > abc. The final result will be expressed in terms 
of matrix elements between state vectors characterized by their particle content and the 
four-momentum each particle has. Such states are improper vectors in an appropriately 
defined Hilbert space, and they owe their name to the fact that they lead in some context 
to undefined mathematical expressions. To avoid problems with improper vectors one can 
consider wave-packets instead of states with definite momentum. A slightly less expensive 
way is by normalizing the vectors with definite momentum with respect to a box of finite 
volume V; and by making explicit that the interaction that causes the decay is not present 
for an infinitely long time but only for a finite time r. 

Both wave-packet and box approach have their respective advantages. In the box nor- 
malization the phase space factors are more clearly accounted for. To use the box normal- 
ization for the decay distribution of an instable particle faces, however, the inconsistency 
that the energy of the particle is assumed to be definite whereas a particle can only decay 
if it is not in an eigenstate of energy. In that respect it is more adequate to represent the 
initial state as a superposition of momentum eigenstates with nearly but not exactly the 
same eigenvalues, i. e. as a wave-packet. 

3.2.1 Volume normalization 

The state vectors \B^,pb) are considered as plane wave functions 

^{t,x) =f^e'^'e-'P^. (3.6) 

The absolute squares of these wave functions represent a probability density and have 
therefore to be normalized to unity. I will do this with respect to a volume V = L^. 
Therefore the method shown here is also known as box normalization. Choosing the volume 
to be a box makes the calculations easier. The idea of normalizing states to a finite volume, 
however, does not only apply to this special choice of volume. The normalization constant 
M is determined by 

1= [ d^x\^{t,x)\'^ =X'^ [ d^x = V (3.7) 

to be 



V Jv 



In a box of finite volume V the three-momenta are discrete if one imposes periodic boundary 
conditions 

^p(t; (0, X2, X3)) = %(t; (L, X2, X3)), (3.9) 
^p{t- (xi, 0, X2)) = ^'p(t; (xi, L, X3)), (3.10) 
^p{t; {xu X2, 0)) = ^p(t; (xi, X2, L)). (3.11) 



3.2 Cross sections and decay rates 
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Then one has, for instance, 



1 ... 1 



V vv 

= e-^^i%(t;(0,X2,X3)), (3.12) 



which requires the eigenvalues for p to be quantized, i. e. labeled by a triplet of integers 

n e Z^, 

p = —j-n. (3.13) 

Because the plane wave functions are normalized to unity, the scalar product between the 
momentum states in the box is a Kronecker delta, 

{p{n')\p{n)) = / d^x^!H'{t,x)^!l{t,x) 
Jv 

= ^ j^d^x exp{-i{p{n') - p{n))x) = S^.^. (3.14) 

Be |/)box and |z)box a state of nif {rrii) particles each with three-momenta pmj {Pm,), 
normalized to the box of volume V. (Multi-particle states are direct products of one- 
particle states.) The transition probability from |z)box to |/)box is given by 

Aw{f, i) = |box(/|5'|z)boxP = |box(/K)box + «box(/|T|z)boxP = | box (/ |T| z)box (3.15) 

The probability for a transition of |i)box into any momentum configuration is 



wif,t) = m2\^oAf\m 

rrif Hmj 



box 



|2 



n ^^|box(/|T|z)box|' 



((^) [n/^V I lbox(/|T|2)boxr, (3.16) 



where in the last step I performed the limit L —>■ oo, using Riemann's definition of the 
integral and equation ()3.13|) . What is the appropriate expression for box(/|^K)box in the 
limit of large volume (i. e. L ^ oo)? We want the continuous momentum eigenstates to be 
normalized as 

(/'!/) = n(2^)'2i^-/5'(p4, -Pn.,) 
m f 

(3.17) 



Jj2E^^. / d^x exp(-z(p'^^ -Pnif] 



X 
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and similarly for In the limit of large volume the right-hand side of equation ()H.17|) 



IS 



X] 



mf 

= V"^ ( n2^-/ box(/'|/)box. (3.18) 



It follows that 



l/>-=lni/#T7|l/>. (3.19) 



and similarly for \i) such that 



box 



. rrii 



and therefore from equation ()H.l(j|l 

T conserves energy and momentum. T is zero whenever the energy and momentum of the 
final and initial state are different. It is therefore suitable to consider the multi-particle 
states 1^) (and |/)) as direct products of eigenstates of the total four-momentum and a 
vector which contains all the remaining characteristics of the multi-particle state, such 
as the masses, three-momenta and quantum numbers of each particle in the state^, 

K) = btot)®IV^.)- (3.22) 

Since T conserves the total four-momentum it acts as the unit operator in the space of 
eigenstates of total four-momentum and thus takes the following form as a direct product, 

T=10M. (3.23) 



^The vectors lip) span the little Hilbert space, see [23 P- lUff-] 



3.2 Cross sections and decay rates 



21 



The T elements between final and initial states then look like 

{f\n) = {£l\l\pSt){i^f\M\ij^) = 6\p,-pf){ijf\M\ij^) ^ 5\p.-pf)Mf,. (3.24) 
In terms of M. we have 

- (rjU) (n / (3.25) 

To make sense of the two delta-functions we write for one of it its integral representation 
in the volume and time interval that is large but explicitly accounted for. To be explicit 
in this subtle point I consider the delta-functions explicitly as distributions in a space of 
test-functions /(p/). 

d'pj[6\p,-pj)ff{Pf) 

= j (fpf 5\pi - Pf)j^ exp(-z(pi - Pf)x)f{pf) 



(27r)4 Jy^ ' (2vr; 



Vt 

d'PfJ^4^'(P^-Pf^f(Pf^- (3-26) 



In the usual short-hand notation this reads 



Vt 

[S\p^-Pf)r = J^,nP^-Pf)■ (3.27) 



Equation ()3.25|) then yields 

2E„, i (27r)3-»-4 



w 



1 \ T/l~™i 

(/' ^) = n ^ T^J^;^. - / dQ l-M/.r, (3.28) 



with 



^g-^%-P/)n^ (3-29) 



the phase space volume for a final state containing nif particles. Usually the quantity of 
interest is the rate, i. e. the transition probability per interaction time r 

TiM-^- (n ^] ^ / (3.30) 
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rrii = 1 (decay) 

In the case of = 1 equation ()3.28|1 reads 



w{f,z) = ^2nr j dQ\Ms,\\ (3.31) 



According to this equation the transition probabihty for decay should increase proportional 
to the interaction time r. However, when the particle has decayed the probability of decay 
should be zero. So equation ()3.31|) is only valid for r much less than the mean lifetime of 
the decaying particle. This conflicts with the assumption of very large r when representing 
the 5^ functions as an integral over space-time. Also formula 1)3.311) contains a factor 2E, 
where E is the energy of the decaying particle. However, as mentioned at the beginning of 
this section, a particle only decays if it has no sharp energy value. These problems can be 
avoided when the initial particle is represented by a wave packet, see section 13.2.21 

rrii = 2 (two-body scattering) 

With rrii = 2 equation ()3.3()|1 is 

The cross section is defined as rate per flux. In the cms (pi = —'P2 = pf)i/\f>i\) the flux is 
given by (p: particle density, v: relative velocity between beam and target particle) 

V \Vi-V2\ \pi/Ei-p2lE2\ p{Ei+E2) 
^ = P''=V = ^:^ = V = E,E2V ' 

and 

s = {pi+ P2f = {El + E2f - {pi + P2f = {El + E2)\ (3.34) 

such that the cross section for scattering of two spinless particles into a final state containing 
mf particles characterized by a set of quantum numbers A reads^ 

^a(2 - rnj) = = _J__ J dQ^^ \m^^^^^\\ (3.35) 

The total cross section is the sum over 

• all numbers of particles M{^/s) that can be produced at a given cms energy ^/s, and 

• for each number of particles all (accessible) set of quantum numbers A, 



^This formula for the cross section is valid in all reference frames where the momenta of the beam and 
target particles are parallel or antiparallel and also in the laboratory frame, see (23 P- 154]. 
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M 



(Ttot = o-(2 anything) = ^ ^o"a(2 nif), 



(3.36) 



mf=2 A 



where M{s) is maximal number of particles that can be produced at a given center of 
mass energy ^/s. Using "anything" in this sense, I adopted the following normalization 
and completeness relation (cf. eq. ()3.17|) ). 



(pi, . . . Amjp'i, • • . = S^fm'^. JJ(27r)^2E„^.(5^(p'„^. - p,r,^)6, 

>=En 



A,n^,A^' (3-37) 



) Pnifj ^nif) \Pli ■ ■ ■ ■> Pnif ; -^r. 



■If THf 



(3.38) 



E 



Ptot,mJPtot)(Ptot| ® / #|^/')(^/'| 



The total cross section then can be written as 



O"tot 



16n 



mf=2 A 



(3.39) 



We further have with pj and pi the total four- momentum of the final and initial state 
respectively (cf. equation ()3.52p and following) 



XI 5Z / ^^™/ = / d'^Pf^^iPi-Pf) ^Q^y 



mt=2 A 



rrif =2 A 



Therefore, 



O"tot 



rf/?|(/3|TK)p. 



(3.40) 



(3.41) 



3.2.2 Wave packets 



To obtain the relation between the respective T-matrix elements and the decay distribution 
of Ti^-K^K^ I follow the lines of j2Zl p. 140ff.], which they use to define the scattering 

cross section of two particles. 

Let la) be the state vector that represents a collection of i?+'s with a four- momentum 
distribution 0(p_b). 



I«) = j d'pBdipl - ml)9{p%)cj>{pB)\B^,PB). 



(3.42) 
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If we consider a collection of fi^'s at rest, (f){pB) will be narrowly centered around pb = 
(m^, 0,0,0), where rriB is the B mass of around 5.279 GeV. The four-momentum eigen- 
vectors \B~^,pb) are normalized as 

{B-^,p'\B+,p) = 2p%6%p-jf), (3.43) 

and the norm of the state \a) is the time derivative of the particle number density 

Nb= / dt{a\a), (3.44) 



where Nb is the number of B^'s produced in the experiment. 
5"-, T- and A^-matrices are defined by 

S=l + zT = l + i5\pB - Ptot)M, (3.45) 

where ptot = EjPi- 
The transition amplitude for the S to decay in a three-body state 1/9) = \7r~ , pi] TT^ , P2; , p^) 
of definite momenta of the three decay products is given by 

{f3\S\a) =i{f3\T\a). (3.46) 

The probability for a decay into a volume of momentum space is the square of this ampli- 
tude integrated with the appropriate measure (see |2Z1 p. 141]) 



P{B]Pi,p2,P3) 



d^pi d^p2 r'^^'' d^p 



X 



j d'pB d'p'M - ml)e{p%)<j>{pB) (3.47) 

x5(p'/-m|)^(<)0*(PB) 
X {B+,p's\T^a){a\T\B+,pB). 

The measure has to do with the normalization of equation ()3.43|1 . To simplify this expres- 
sion by an approximation it is useful to express it in term of elements of Ai. We have 

(/3|r|a) = j d^pB5{pl-ml)d{p%)ct>{pB)5\pB-Ptot)mM\B^,PB). (3.48) 

If the momentum distribution (j){pB) is sufficiently narrow, we can approximate the above 
expression by assuming that the M. element is independent of and that it can therefore 
be pulled outside the integral. For the decay probability we then obtain 



P{B;pi,P2,P3) 



d^Pi ^2+^2 d^p2 d^p 



Pi '^Pl Jp2 2p2 Jpi 2^3 

(3.49) 



X mM\B+,pB)? j d^pBd^p'^-mlWiplHiPB 
X Sip's' -ml)e{py)<P*{p'j,)S\pB-ptot)S\p's -Ptot). 
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For further simplification 1 introduce the Fourier-transforms of the momentum distribu- 
tions, 

= j d'pBSipl - m|)^(p°j)0(pB)e*"^^ (3.50) 
After some algebra one arrives then at 



(Ppi ^^+'^2 d'^p2 r^^^' d^p 



'3 



The integrations in the space of final state momenta can be written as 



3 fPi+k. 



i=l 

with 



3 rPt+ki li 

= n / ^^'(Pi +P2 +P3 -Ptot). (3.53) 

i=l * 

This can be checked using a test functions 

/(X^^^'X^^»'^l'^2,P3)) =f{Pl,P2,P3) (3.54) 

i i 

as follows (ptot = Y^^i=iPi, Ei = ^Jnlf + pf)■. 

j d%ot dQsf (Ptot , Ptot , Pi , P2 , P3 ) = 



Pi+ki 



d?p, 



n/. 



=1 ■'P 



P2,P3j 



(3.51) 



n / I d'ptotdQs (3.52) 



3 rpi+ki Ji 

n/ ^M^P2,P3)- (3.55) 



=1 ■'P 

If we integrate over all values of ptot the above equality is not exact. However, we will 
use this equality to evaluate expressions such as that in equation ()3.49|1 . There we assume 
the momentum distribution function 4>{pb) to be narrowly centered around a certain value 
((772^,0,0,0) in the B rest system). Because of the factors S^{pb — Ptot)'5'^(PB — Ptot) the 
integrand in equation ()3.49p vanishes unless ptot = Pi + P2 + P3 takes on values near the 
central value of p^. So the mistake we make by integrating over all values of ptot can be 
made arbitrarily small by making the momentum distributions arbitrarily narrow. 
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We can now make the substitution of equation ()3.55|) in equation ()3.51|) and perform 
the integration over ptot- This gives 



(3.56) 



By Fourier-transforming back $(x), using 

c/V(p' - m2)^(/) = (fpdp^i^p""^ + m")) 



d^pdp^ 



X 



5 (p^ + ^/p^ + mA +6(p°- ^/p^ + mA 9{p' 



(3.57) 



d^p 
~ V 

and assuming once more that the momentum distributions are narrow one can show that 
where the last equality is equivalent to equation ()3.44p . The particle density p{x) is then 



p{x) 



pI 



such that 



B 



d'^x p{x) 



(3.59) 



(3.60) 



is the total number of B's produced in the experiment. 
Thus equation ()3.56|) reads 



P{B;pi,p2,P3) 



B 



2TXp 



B 



dQ3mM\B+,PB) 



(3.6i; 



The differential number of decay events that yield decay products in the momentum volume 
dQ^ can be read off to be 



dP{B;pi,p2,P3) Nb 



dQ3 



27rp 



, mM\B\pB)\'. 



(3.62) 



In this way the A^-elements between momentum eigenstates determine the decay distri- 
bution of the B meson that is characterized by <f){pB) as a wave packet. 



3.3 Three-body phase space in a Dalitz plot 



27 



3.3 Three-body phase space in a DaHtz plot 

In this section I try to show, that the phase space volume J dQ^ is constant over the Dahtz 
plot. There are alternative proofs in the literature, e.g. [3 p. 140f.] and |27| p. 159ff.]. 
However, to me there are some rather involved steps necessary. Therefore I try here my 
own derivation which is still similar to the one of ref . j2Z| . 

The factors J can be related to an integration over the four-momenta with a 6- 
and a 6 function that restrict the integration to the respective mass shells and to positive 
energy values {f{pi): test function, cf. [SI p. 67] and j221 p. 497]): 



d%5{pi-m1)e{p^)f{p, 



d'pidp^Sipf - (p1 + m2))0(pO)/(pO,p-; 
d'^Pi dp^i 



X 



m- 



+ S[p^ 

d^Pi 



p^ + mf 



(3.63) 



1 



d^Pi 
2E, 



2^/pf + ml 
f{Ei,pi), 



fi^fi+mlp) 



where Ei = \/^~+ 

A decay of type -k^-k^K^ is completely characterized by the four particle masses, 

the three-momenta and the energies. All these quantities can be expressed in just two scalar 
variables. In chapter |3] I did this explicitly with the variables Si2 and S13. Since these two 
variables are related to and E2 by 



s,, = ml + ml-2msE, and 



Sl3 



m 



B 



TTln 



2mBE2 



(3.64) 



the latter suit equally well. Whether one chooses either pair to define the axis of the Dalitz 
plot is mere convention. To show that phase space is constant over the Dalitz plot it is, 
however, more convenient to work with E2 and E3. 

Since the decay is completely specified by just two variables, also the matrix elements 
of can be expressed in terms of these. 



{P\M\B^,Pb)=M{E2,Es). 



(3.65) 



In the expression for the number of events in a certain region of phase space, equa- 
tion 1)3.611) . we can therefore pull ^^(£'2, -Es) outside the integration over the remaining 
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three variables dq, with dQ^ = dE2 dE^ dq. So let us evaluate the integrals over the re- 
maining variables. I will do this with respect to the rest frame of the B. We start with 
(see equation ()3.53|) ) 

JdQ3=flJ ^f{pi+P2+P3-PB). (3.66) 



i=l 

The factor we replace with 

d%Sipl-mj)9{pl), (3.67) 
see equation ()3.63|1 . For the factors and we make the substitution 

^^ |KM(|P3lV(cos.)# ^ i|fi|*os»)*<;E3. (3.68) 

(3.69) 

where without loss of generality we have chosen a coordinate system such that the polar 
angle 6 coincides with the angle between p'2 and P3. This will be of use later. Equation ()3.66|) 
now reads 



JdQ3= J d^pi6{pl-ml)e{pl)S\pi+p2+P3-pB) 



X ^d(cos^) dcj) \p3\dE3 ^d{cos9') dtp' \p2\dE2. (3.70) 

The integration over 0, 0' and 9' yield a factor Svr^ since the integrand does not depend on 
either of these variables. Performing then the pi integration and the trivial integration 
yields 

j dQ-s = 271^ j 5[ml + {E2 + ^3)' - 2mB{E2 + ^3) - \P2? - \pz? - '2\p2\m cos9- ml] 

X 9{mB -E2- Es)\p2\\p3\d{cos9) dE2dEs. (3.71) 

Now the integral over the variables other than E2 and £'3, from which only M depends, is 
reduced to an integral over cos 9, 

j ^6[ml + {E2 + ^3)' - 2mBiE2 + E^) - Iplp - jp-gp - 2\p2\\n\ cos9 - mj] 

X \p2\\p3\d{cos9) 

1 

5[cos^-2o]c/(cos^), (3.72) 



2 J-i 



3.3 Three-body phase space in a Dalitz plot 
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with 



ml + {E2 + E-iY - 2mB 


{E2 + ^3) - \P2 


' - \P3 


— m\ 


2 


P2I 


|P3 





(3.73) 



The integral then gives 



^^(1 - ^0) = ( 1 



m| + (E2 + E3)'-2mB(^2 + 


^3) - \P2 




— m\ 


2bl||p3 





7^0 [2IP2IIP3I -ml- {E2 + ^3)' + 2mB(^2 + ^3) + \P2? + \n? + m?] 



2\ {El - ml){El - mf) + ml - ml - ml - ml - 2E2E3 + 2mB{E2 + ^3) 



(3.74) 



Altogether, we have found that 
J dQs = 7r^ J e{mB-E2-Es 



2J{El-ml){El-ml 



m 



B 



m2 — m^ 



2E2E3 + 2mB{E2 + E3 



dE2 dE-i. 
(3.75) 



The theta functions define a boundary in the E2-E3 plane outside which the phase space is 
zero and inside which it is constant vr^. The first theta function gives a straight line in the 
E2-E3 plane. The more interesting characteristics of the Dalitz plot boundaries are given 
in the second theta function. The condition that its argument is zero corresponds to the 
energy-momentum conservation in the special configuration of cos^ = 1, as can be seen as 
follows.^ 



Pi +P2 +P3 = 0, 

Ei + E2 + Es= mB. 
The first equation is more useful in the form (cos 6' = 1) 



P2^ +P3^ + 2\p2\\p3\ 



Pi 



(3.76) 
(3.77) 



(3.78) 



Expressing all momenta in terms of E2 and E3 and the masses, using the energy conserva- 
tion condition yields 



El 



mi + Ei 



mn 



+ 2J{El-ml){El 



{mB - E2- E3 



m 



l5 



(3.79) 



which is equivalent to 



2J{El-ml){El 



mn 



m^ 



mn 



ml + 2mB{E2 + ^3) - 2^2^3 = (3.80) 



■^Also cos 6* = —1 corresponds to events on the boundary of the Dahtz plot, see figures 
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where the left-hand side is the argument of the theta function in equation ()H.75|) . The 
points where this equation is satisfied he on the closed curve shown in figure IH.ll 

That the configuration with cos 6 = 1 gives the boundary of the Dalitz plot was reflected 
in the above calculation in the transition from the delta to the theta function, i. e. from 
equation ()3.72|1 to equation ()3.74|1 . 



3.4 Elastic scattering of two spinless particles 

Scattering of two particles, say a and b, where the final state consists also of two particles a 
and b is called elastic (two-body) scattering. In this case the A^-matrix elements depend at 
most on the two independent scalar variables s (total energy squared in the center-of-mass 
system (cms)) and t (4- momentum transfer). 



S = (pi +P2f 

t={Pi- Ps) 



(3.81) 



where the particle a (6) is labeled by "1" ("2") before the scattering and by "3" ("4") after 
it, and the squaring is meant to be with respect to the Lorentz-invariant scalar product 
p2 — (j)0)2 _ p2 _ (^p0)2 _ (^pi)2 _ (^p2^2 _ (^p3)2_ ^]^g matrix element does not depend 

on more than these two kinematic variables has the following reason: The matrix element 
is a Lorentz-invariant quantity, i. e. a scalar. It can therefore be only a function of scalars. 
The only scalars that can be formed out of pi and p2 are s and t j2Zl P- 168]. One can 
convince oneself also by expressing explicitly any kinematical quantities in terms of these 
two, similarly to what I have done for the kinematics of a three-body decay in chapter |3] 
The momentum transfer t and the scattering angle 6 (in the cms) are related (m^ = 
rrih = m, q: cms momentum): 

t = 2ml-2El+p,ps 
= 2j92(cose-l). ^'^■^ ' 

The A^-matrix elements may thus be written as a function of s and z = cos 6, 

{ab\M\ab) = M{s,z). (3.83) 

The differential cross section for elastic scattering, da{ab ab), per solid angle dQ = 
dzd(j) is given in terms of the A^-matrix element. 

It is sometimes useful to work with the amplitude 

f{s,z) = -^M{s,z), (3.86) 
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Figure 3.1: Boundary of the Dalitz plot according to equation ()3.80|) . with the labehng 
1 ^ 7r~, 2 ^ 7r+ and 3 ^ . The s's are given by equation ()3.fi4j) . 
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such that 

(3.87) 

3.5 Unitarity 

Collecting all quantities characterizing a free-particle state in one single symbol a (or /3 or 
the like), cf. equation ()3.37p . we have the following short-hand notation for the normaliza- 
tion, 

ip\(5) = b{a-(5). (3.88) 
Using the completeness relation ()3.37|) we find 



b{a - 7) = out(tt|7)out = y C?/3out(a|/5)inm(/9|7)out 

= y" diiS^fiS\^ = \SS\^. (3.89) 

Applying in a similar vein the completeness relation using the out-states one obtains 

[5'^^]a7 = 5(«-7). (3.90) 

These two relations mean that S is an unitary matrix. 

In section 13.21 1 showed in some detail the two methods to deal with state vectors 
normalized to delta functions. I will often assume that one can anyhow discretize the states 
and normalize them to unity for instance as done in section l3.2.1l Then, the probability that 
a system in an initial state \i) evolves in any final state is one. This postulate is sometimes 
called probability conservation. The probability for the transition |a)in — > |a)out is given 
by the absolute value squared of the corresponding S-matrix element. The condition for 
probability conservation thus reads 

E = = {^\S^S\z) = 1. (3.91) 

/ / 

This means that the diagonal elements of the operator S'^S are 1. The off-diagonal elements 
are zero as can be shown by the following argument. 

Take two orthonormal state vectors \m) and |n), and set the initial state \i) to (|a;p + 
|/?P)~^/^(a|m) + P\n)) with arbitrary a and /5. Eq. ()3.9H) then implies 

'amm\S^S\m) + \P\^{n\S^S\n) 



+ a*(3{m\S^S\n) + S^S\n) + a(3* {n\S^ S\m)) = 1. (3.92) 
Since according to eq. (j3.91|l the diagonal elements are 1 we get 

a*(]{m\S^S\n) + ap* {n\S^ S\m) = 0, (3.93) 
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and because a and /3 are arbitrary it must be that 

{m\S^S\n) = {n\S^S\m) = 0. (3.94) 

All in all we have derived that S'^S is the unit operator, 

S^S = 1. (3.95) 

Multiplying eq. ()3.95|) from the left with S and from the right with we get that also 
SS"^ is the unit operator. This means that S* is a unitary operator. 
Plugging S = 1 + iT into eq. (|3.95j) yields 

T^T = i{T^ - T). (3.96) 



3.6 Partial wave amplitudes (without spin) 

The mathematical objects that are interpreted as elementary particles are irreducible rep- 
resentations of the Lorentz group. For the question whether in scattering of two particles 
resonances or new particles are formed by these, it is therefore important to decompose the 
two-body states \ip), which are in general not an irreducible representation of the rotation 
group, into states of an irreducible representation. In other words we have to expand it in 
eigenstates of angular momentum. 

\i^) = Y,CjM\JM,k) (3.97) 

Such an expansion is also useful inasmuch as angular momentum is a conserved quantity, 
such that the matrix element can only be non-zero between states of the same total angular 
momentum and third component, i.e. we can write 

(J'M',A'|A^|JM, A) = 5rj5M'MMi,^. (3.98) 

As a consequence of the Wigner-Eckart theorem A^a'a does, as the notation suggests, not 
depend on the third component M of angular momentum. If we expand in- and out-states 
in eigenstates of angular momentum we thus obtain^ 



{^l;J\M\i^^)= Yl CjMC},,,,{J'M'J\T\JM,t) 

J,M,J',M' 



J,M,J',M' 



(3.99) 



J,M 



l-Y,{'^J+l)Pj{z)Mi^, 



Air 
J 



*For details see e.g. [23 p- 177fF.]. Pj: Lcgendre polynomials. 
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Note that {tp f\Ai\tpi) is a function of s and z = cos6, while A^a/a is only a function of s, 

^) = + ^)PA^)M'{s), (3.100) 

J 

i. e. the 6 dependence is entirely contained in the Legendre polynomials Pj, and the depen- 
dence on the cms energy is contained in the partial wave Al-matrix element Ai'^{s). Again 
(cf. equation ()3.8(ij) . we may define a sometimes more convenient quantity, the partial wave 
amplitude 

such that^ 

f{s, z) = 5^(2J + l)Pj{z)f{s). (3.102) 
J 

The inverse of equation ()3.102|) is 

f'{s) = ^l\j{z)f{s,z)dz. (3.103) 

3.7 Phase shifts 

The postulate of probability conservation says that an initial state is to be found with 
certainty in any of the final states. By definition of the concept of a conserved quantity, 
the probability of finding the system in a final state that differs in a conserved quantity 
is zero. Two conserved quantities that characterize hadronic states are the total angular 
momentum and isospin.^ In hadron spectroscopy it is therefore adequate to decompose in- 
and out-states in eigenstates of these two conserved quantities. As a consequence of the 
conservation laws unitarity relations hold for each component. 

As in the discussion of the partial wave expansion we can define a new scattering matrix 
element by factoring out Kronecker deltas representing the conservation of total angular 
momentum and isospin and their respective third components (cf. eq. fj3.98|) ). 

{J'M'rr,,A'\M\JMIh,A) = 6j,j6M'MSi'i6j-j,Mi (3.104) 

where A symbolizes all quantum numbers other than angular momentum and isospin. A4j 
depends in general on A; I dropped the indices for brevity. As consequence of the Wigner- 
Eckart theorem it does not depend on the third components of angular momentum and 
isospin. 

It is possible (cf. jSl p. 152]) and useful to define an S operator in the little Hilbert 
space where sofar we only have defined the elements of Ai, such that 

{J'M'I'I'„ A'\S\JMIh, A) = 6j,j6M'MSri6r^j,Sf, (3.105) 

have dropped the indices "A" for brevity. 
^Isospin is only approximately conserved in the strong interactions. 
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and 

^/(A'A) = 5j,jSM'M5ri5i'^i,5A,A - zA^/(A'A). (3.106) 
The condition of probability conservation (eq. (j3.9H) ) then takes the form 

1= J2 \{J'M'I'I'^,A'\S\JMIh,A)\^ = \Sj^{A'A)\\ (3.107) 

From this one can derive in a similar vein as the unitarity of S in general has been derived 
(section l!?3|) that also each partial wave ^'-operator, Sj, is unitary and that J^j satisfies 
the same constraint as Ai, 

(SfySf = 1, (3.108) 
{MjyM-l = i{{M-{y - M^j). (3.109) 

If only the elastic channel is open, i. e. A' = A, then S*/ is just a complex number. Unitarity 
is then the condition that this complex number be of modulus unity and so can be written 
in the general form 

Sf = e^'^', (3.110) 

where we have introduced the phase shift 5 G M. (The factor 2 is convention.) 

If also other than the elastic channels are allowed, it is useful to introduce another 
(real) quantity that characterizes a certain scattering process, the elasticity rjj. This is in 
a given partial wave the square root of the fraction of the probability of elastic scattering 
and the total probability of any scattering, which is unity. 

With this definition, the elastic partial wave amplitude S'/(AA) can be written in terms of 
the elasticity factor and the phase shift 5j in the respective partial wave characterized by 
angular momentum and isospin, 

Si{AA) = 7]ie^'^i . (3.112) 

In the case of purely elastic scattering, i. e. when no other final channels are available as 
the initial one, we have rjj = 1. Using this form of Sj{AA) the partial wave A^-matrix 
element for elastic scattering takes the form 



For elastic scattering we have {rjj = 1) 

A^/(AA) = 2e'^' sin 5/. (3.114) 
If the elasticity tends to zero, rjj = 0, the corresponding partial wave A^-element tends to 

Mi{AA)=i. (3.115) 
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3.8 Resonances 

3.8.1 Argand diagram 

The relation (see eq. (jH.llHj) ) 



Mj{AA) = Mi = (3.116) 



can be transformed in 



Mj = i- ini cos 25/ + r// sin 25/. (3.117) 



From this we can see that the complex number M.j can be represented as a point in the 
complex plane in the interior or on the boundary of the unitarity circle. 

3.8.2 Breit-Wigner resonance 

The epitome of a resonance is a peak in a partial cross section with definite total spin; a 
peak that is caused by an amplitude with elasticity 77"^ = 1 and a rapid shift of the phase 
with energy through 5'^ = tx /2. In the Argand diagram this is represented by a rapid 
movement through the top of the unitarity circle, see figure ESI (a). 

Such a resonance is called Breit- Wigner after the inventors of the formula for the cross 
section describing such a resonance. For the scattering amplitude we have in the Breit- 
Wigner case (dropping isospin labels) in the neighborhood of m|j, the mass of the resonance 
squared, 

with r the width of the resonance. The phase shift is then given by 

r/2 

taii5=^ . (3.119) 

mjj — s 

If s changes from values below m\ to values above, 6 changes rapidly from zero through 
7r/2 to vr (or odd multiples thereof). 



3.8.3 Background 

In general resonances are as in the Breit-Wigner case characterized by the two parameters 
r and m/j and a rapid phase shift. However, these quantities are in general not directly 
related to the position and width of a peak in cross section. They rather denote the 
position of the pole oi M..^ Also the rapid phase shift does not necessarily occur at 7r/2 
(or a multiple thereof). It may happen that the rapid phase shift begins only when the 
phase has already appreciable values different from zero. 



3.8 Resonances 



37 




Figure 3.2: Argand diagram, 
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(a) Breit-Wigner; the amplitude is maximal (b) Resonance with background phase of 

at s = m|j. about 7r/2. The amplitude is zero at s = 

771%. 



Figure 3.3: Resonances in an Argand diagram. 








FIGURE 13.3. Four possible resonances. The 6(p) plots show the resonant phase 
shifts for djjg = 0, 77/4, n/2, and 3nl4. The sin^ 6(p) plots show the corresponding 
behavior of the partial cross section [apart from a factor 47t(2/ + l)//>^]. 



Figure 3.4: Phase shifts and corresponding characteristic feature of the cross section; from 
p. 242]. 



3.9 Optical theorem 



39 



Figure EiU shows four possibilities of a resonant phase shift and its corresponding char- 
acteristic behavior of the cross section. Appreciable background phases are not just a 
mathematical possibility but seem to occur in actual physical reactions. One prominent 
example is the /o(980) which appears as a dip in the tttt spectrum because the correspond- 
ing phase shift increases rapidly through vr and not through 7r/2 as for a pure resonance. In 
the Argand diagram the amplitude moves rapidly through the bottom of the unitarity cir- 
cle, see figure Eini (b). The experimental data and a parametrization of the relevant phase 
shift is shown in figure 13.51 together with the inelasticity. In figure 13.61 the cross-section of 
the isoscalar S wave is compiled. Such a structure can be explained by two rapid phase 
shifts through vr and 27r, as I roughly simulated in figure l37fl 



3.9 Optical theorem 

In the special case where the initial an the final state are the same, so-called forward 
scattering, the unitarity condition for the T matrix (equation ()3.96|l ) reads 



2ImT„„ = j dl3\Tp^\\ ^^-^^O) 
ImTaa = Svr^pv^cTtot- 

To obtain a form of the optical theorem in terms of }A we have to take care again 
about delta functions. Therefore, we cannot substitute simply T = S'^{pi — Pf)M. in the 
above equation. In ref. jHl p. 147] this is done by dividing through delta function, which 
I want to avoid. The alternative that I prefer is again by means of wave-packets, cf. ^TJ\ 
p. 183ff]. I consider states of the form 



= d''pM'^{p^)\p^) ® I^A.). (3.12i; 



Sandwiched between two such states \i) and |/) the operator equation (j3.96p takes the 
form 



/ d^pM'^{p^)^^^^*{p^){Mf,-M], 



M „ 

d^p,$«(p,)$(^)*(p,) dQ.^^ Ml^^^^M^^. (3.122) 

mf=2 A 

Since the momentum distribution functions $ are arbitrary the integrands have to be equal. 
Therefore, by equation ()3.39j) . we obtain for forward scattering, i.e. 2 = /, 

ImMii = STT^pv^atot, (3.123) 



^For a discussion of pole parameters and Breit-Wigner parameters see |2()| . 




Figure 3.5: Resonant phase shift starting from an appreciable background phase in the 
isoscalar S wave. In the lower panel the behavior of the inelasticity in this partial wave is 
shown; from |29j . 
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Figure 3.6: Compilation of the isoscalar S wave. Three major peaks or alternatively two 
major dips can be seen. Copied from jT2j, see references therein. 
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Figure 3.7: (a) Very rough parametrization of a phase shift that gives rise to a cross-section 
comparable to that of figure ESI (b) sin^ of the rough parametrization of the phase shift, 
which gives the characteristic form of the corresponding partial cross section. 



3.10 Diffraction peak 
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the right-hand side of which is the same as in equation ()H.12()|) . in this sense the delta- 
function relating T and ^A has canceled. The result of equation ()H.12Hj) could have been 
found also by first showing that M. satisfies the same form of unitarity condition as T, i. e. 

i{M^ -M)=M^M, (3.124) 

see, for instance [3 p. 147]. 

Because Ai and T satisfy the same form of unitarity condition (equations (j3.9fi|l and 
fl3.124j) ) I will often use T and Ai interchangeable. In particular when discussing unitarity 
constraints (section 15.2.11) I will use T elements without worrying much about the relation 
to the Ai elements in my conventions. None of the results obtained later hinge on the 
difference between T and TW. Thus I adopt the notation of ^JEDj for formulating unitarity 
constraints. 



3.10 Diffraction peak 

The element of phase space of a two-body, i. e. m/ = 2, final state is 

dQ2 = ^^^'(Pi +P2-Pa- Pbl (3.125) 

where pi and p2 are the four-momenta of the two final state particles, and Pa and pb the 
four-momenta of the two initial particles. We now discuss this expression in the center-of- 
mass frame where we have, by definition, that Pa + Pb = 0. Thus, cf. [HI p. 139ff.], 

dQ2 = ^^^'(Pi + P^MEi + E,-E,- E,) 

1 <Pp2 



I 2E2 

2ymf+p2 
1 \p2\^ d\p2\d(l) dd sin6' 



5Uml+p2' + E2-Ea-Eb 



'2\jm{^P2 2^m^ + p2' 
For the delta function we can use 



b{^lm\\p2'\^lml\p2'-E,-Eb). (3.126) 



with q the zero of /(IP2I) and find 

5(V'^?+P2'+ 4^ni^^ -E,- E,) = 5(|p2| - q)j^^^, (3.128) 

with 

^[{Eg + Eb)2 - mf - mi]2 - 4mfmi 
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Thus we obtain the relation between the two-body phase space and the scattering angles 
9 and (f) in the cms frame: 

We can now make an estimation that reveals a certain behavior of the total cross section 
for two-particle to two-body scattering, cf. equation ()3.35p . We assume that the matrix 
elements of M. are smooth functions of the scattering solid angle Vt and that they fall off 
to zero for large angles. Then we can find a solid angle /STL such that 

dn\Mu^)i\^ >]^m\Mii\\ (3.131) 



Then we obtain 



atot > ^a(2 ^ 2) = Y^J^ j dQm.^ 



1 q 



dftlM 



{/A)i| 



|2 



b) 



- 167r2pv^4(E, + Eb) ' 

1 q 



Afi647rVsc^tof (3.132) 



IGtt^P^ 4{Ea + Ek 
In the cms p = q, Ea + E^ = a/s. So we obtain an upper bound for Afl, 

An < —\ . (3.133) 

Hadronic cross sections show a diffraction peak, i. e. the differential cross section falls 
off exponentially with t = 2p'^{cos9 — 1) (cf. eq. ()3.82|1 ). In the special case of a totally 
absorbing disk of radius R (see e. g. jHH p. 136f.] and |2Z1 p. 210f.]) we have approximately 

oc exp ^ . (3.134) 



dt ' \ 4 ^ 

In figure 13.81 the differential cross section for elastic proton-proton scattering is shown. 
There we see that already for incident momenta of 5 GeV (i.e. ^/s ~ 3.4 GeV) the fall-off 
is visible but that the diffraction peak is more and more pronounced for higher incident 
momenta (in the figure up to 1480 GeV, i. e. y/s ^ 53 GeV). 

As shown in figure 13. 9[ typical hadronic cross sections rise with energies for energies 
greater than roughly 20 GeV. Thus also equation ()3.133j) shows that the forward peak in 
the two-body to two-body total hadronic cross section becomes more and more narrow at 
these high energies. 



3.10 Diffraction peak 



45 




Figure 3.8: Differential cross section for elastic pp scattering as a function of the square of 
the momentum transfer, \t\] from jHII p. 137]. 
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Figure 3.9: High energy behavior of typical hadronic cross sections; from [T^ . 



Chapter 4 

Kinematics of a three body decay 



In this chapter I collect kinematical relations that could eventually be used for a partial 
wave analysis. As a numeric example I take the decay ti~tt^K^ . 



4.1 Rest system of decaying particle 

First I consider the kinematical relations that hold in the rest system of the B meson. By 
definition its three momentum is zero, ps = 0, in this frame of reference. Its energy is 
given by its mass triB- The squares of the respective two-body invariant masses are defined 
by 

S7r-7r+ = (Ptt- +V^+f = {En- + E^+f " (Pn- + Pn+f , 

Sn''K+ = (Ptt- +PK+f = {E^- + Ek+Y - {Pn- +PK+f, (4.1) 
Sn+K+ = (P7r+ +PK+Y = {En+ + Ek+Y - {p^+ +PK+Y- 

Recall that these quantities are Lorentz invariant; the energies and the three-momenta for 
one particle pair can be taken with respect to any reference frame. (But it has to be the 
same reference frame for the energy and the three- momentum.) For the sake of brevity 
I will sometimes refer to these quantities just as "two-body masses", "pair masses" or 
the like. These masses (squared) of one pair of particles are related to the energy of the 
remaining particle, e. g. 

= (niB - Ek+Y - (4.2) 
= m| + m|-+ - 2mBEK+, 
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and therefore also to the modulus of the three-momentum of the remaining particle: 



-e 7-12 2 
Pk+ = Ek+ - m 



K+ 



m% + m|.+ - Sr, 



2ms J ""^^ (4.3) 

m% + + sl_^+ - 2{m\m\+ + m%s^-^+ + m^+s^-^+) 



4m^ 



Similarly, one obtains expressions for E^^- and and Ej^+ and in terms of s,r+i<'+ 
and respectively, see table 14. ll 



4.2 Two-body system 

We now make a Lorentz transformation into the center-of-mass system (cms) of the two 
pions, say. The kinematic situation in this reference system is shown in figure EHl Compare 
it to the kinematic situation in scattering in the cms shown in figure I1I21 The appropriate 
transformation from the B rest system to the tttt system is a Lorentz boost along the 
direction opposite to Pk+- The relative velocity between the B rest system and the two- 
body system considered now is (in unities of c and with respect to the direction oi px+) 



\Pk+\ 



(4.4) 

Ejr- + E^+ ms — Ek+ 
This velocity defines the hyperbolic angle a by which a Lorentz boost along a given axis 
(here the axis given by Pk+) can be specified: tanha = v. The matrix that represents the 
Lorentz boost then reads 

/ cosh a — sinh a 0\ 
sinh a cosh a 





V 





1 
1/ 



(4.5) 



Accordingly, the energy of the B in the new reference frame is 

ms - Ek+ 



E'n = ms cosh a 



ms- 



rriB 



m 



rriB 



B 



m 



K+ 



2m B 



(4.6) 



m 



B 



m 



K+ 



In the center-of-mass system of the two pions we have E'^ 
of the in this system is therefore given by 



E' 



— m 



K+ 



^s^-^+ + E'j^+. The energy 



(4.7) 



4.2 Two-body system 
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B+ 

-m 




Figure 4.1: Kinematic situation for the three-body decay tt^tt^K'^ in the vrvr rest 

system with the two-body angle 6 as defined in equation ()4.13j) . z = cos 6. Top of the 
figure: before the decay of the B. Bottom of the figure: after the decay and the strong 
final state interactions. 



• * * • 

^_ 




Figure 4.2: Kinematic situation for scattering in the center-of-mass system. Top of the 
figure: before the strong interactions. Bottom of the figure: after the strong interactions. 
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By definition of this reference frame the three-momenta of the two pions are equal, = 

So are, by momentum conservation, the three-momenta of the B and the K, ff^ = 
I denote the respective moduh of the three-momenta of the pions and the B (or K) as k 
and /. I is determined by energy conservation in the new reference frame, 

y/m| + /2 = ^mK+ + /2 + (4.8) 
which equation can be solved for /, yielding 

~ As 1 ■ ^ ^ 

As to the three-momentum k of the pions, it can be obtained by evaluating the Lorentz 
invariant quantity St^-t^+ in the new reference system, i. e. 



^/i— 7 = ^ml^ + F + ^ml^ + P, (4.10) 

which upon solving yields 

4s - + • V • / 

The respective energies of the pions in the primed reference system are then 



K- = Vm2_+P, El^ = ^ml^ + k^. (4.12) 

One relevant quantity for a partial-wave analysis is the angle 9 between the respective 
directions of the tt" and the in the new reference frame, i. e. the center-of-mass system 
of the two pions. We are now going to derive an expression for the cosine of this angle, 
i.e. z = cos6, in terms of the two variables of the Dalitz plot, 3^^-.,^+ and Stt-x+j and the 
masses of the four particles concerned. By definition we have 

Suitable expressions for E'^_, E'j^+, k and / are given, respectively, in the equations ()4.12j) . 
dOI), (jOTIl and (gH). For p^-pK+ we have 

Pn-PK+ = K-Ek+ - Pn~PK+- (4-14) 

For Ex+ we have already derived a relation between it and the Dalitz plot variable s^-^+ 
(equation ()4.2j) ). The scalar product between the two three- momenta can be expressed 
by the moduli squared of the three-momenta of all three decay products, by transforming 



4.2 Two-body system 
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the relation for momentum conservation first by putting it in an appropriate form then 
squaring it and finally solving it for Ptt-Pk+- 

P-K- + Pk+ + Pn+ =Pb = 0, 

Pit- +PK+ = -P7r+, 

pt- + p1+ + 2p^-pK^ = pt+ , (4- 15) 



P-n~PK+ 



->1 -Q. ->1 



The squares of three momenta in terms of the three two-body masses are given by equa- 
tion ()4.3j) and the respective entries in table 14. l1 To express all quantities in the two Dalitz 
plot variables s.^-.,j+ and Stj^k+, it remains to find an expression for s.^+k+ in terms of 
these. Such an expression can be obtained as follows. 

S^-^+ + S^-K+ + S^+K+ = {Ptt- +P7T+f + {Ptt- +PK+f + (P7r+ + PK+f 

2,2,2 

= V^- + + Vk+ 
+ pI- + pI+ + p\+ + 2{p^-p^+ + p^-Pk+ + Pn+PK+)- (4.16) 

Because of energy-momentum conservation we have 



fnB=PB = (Ptt- + Pn+ + Pk+)^ 

= pI- + pI+ + Pk+ + 2{p^-p^+ + p^-Pk+ + Pn+PK+)- 



(4.17) 



Using this result and since = etc., equation ()4.1fij) reads 

S7r-7r+ + S^-K+ + S^+K+ = + + + m%. (4-18) 

So to obtain expressions in the two Dalitz plot variables s^-7r+ and S7r-x+ one can replace 
every occurrence of s^+k+ by m| + m^. + m^^ + - s^-^+ - s^^k+. 

The corresponding formulae for the other two two-body systems, i. e. the cms of the 
negative or positive pion and the kaon, can be obtained by making appropriate substitu- 
tions. In the case considered so far, i. e. the cms of the two pions, the kaon played the 
role of the spectator, the negative pion the role of the particle with respect to which the 
angle 6 is measured. In the cms of the negative pion and the kaon, the positive pion is 
the spectator, the reference particle for the angle I take to be the negative pion. (It could 
equally well be the kaon, this is just my convention. With the kaon as reference for the 
angle, z would change sign.) In the cms of the positive pion and the kaon, the negative 
pion is the spectator and the kaon (or the positive pion) is the reference particle for the 
angle. 
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4.3 Pair masses and two-body angles 

4.3.1 Resonance band or angular peak? 

In a typical Dalitz plot for -k^ti^K^ each decay event is assigned a coordinate 

pair (s7r-7r+; With the kinematic constraints discussed in this chapter the decay 

events can be identified^ equivalently by other pairs of kinematic quantities. For the task of 
identifying resonance masses and spins one is interested in the invariant pair masses St^~t^+, 
S-K-K+ and SK+-K+, and in the two-body angles 6*7^-^+, 6.„-k+ and 6k+-k+- In figures l01 to l4.8l 
lines of constant s and z = cos6 are drawn using the Maple command implicitplot. With 
any two out of these six patterns one can construct one of (2) = 15 coordinate systems for 
the three-body decay events. 

The possibility of using coordinate systems consisting of angles and pair masses shows 
that the interpretation of structures in the Dalitz plot event distribution like bands (see 
figures im and is a priori not unique. A concentration of events in a certain region of 
the Dalitz plot can obtain because certain pair masses are more likely to be produced in 
the decay or because certain angular distributions in a two-body subsystem are preferred. 

4.3.2 Partial waves 

The coordinate lines of the two-body angle can eventually be used for a partial wave 
analysis. In the decay -k^tx^K^ the line of constant Ztj-k+ = could be of particular 

importance. In [Hj the reference amplitude and phase is taken to be the signal from 

fi+ ^ if*(892)7r+ ^ K+n-TT+. (4.19) 

K*{892) has spin J = 1 and should thus show its particular angular dependence given by 
the Legendre polynom 

Pl{Zn-K+) = z^-K+ = cos9^-K+- (4.20) 

The partial wave amplitude representing a decay via the i^r*(892) is therefore expected to 
vanish where Zt^-k+ = 0. 

In figure lOl the three lines Zj,-j,+ = 0, Zt,-k+ = 0, Zt,+k+ = are plotted. 



^Actually, the decay events are not uniquely determined by pairs of two appropriate kinematic quantities 
like (s^-7[.+ , s^-j^+). The mirror image of a kinematic situation in a three-body decay is not distinguished 
from its original by such coordinates. 
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Figure 4.3: Lines where 5^-,^+ = (5, 10, 15, 20, 25) GeV^. Abscissa: 5^-,^+ [GeV^]. Ordi- 
nate: St,-k+ [GeV^]. 
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Figure 4.4: Lines where 5^-;^+ = (5, 10, 15, 20, 25) GeV^. Abscissa: s„-^+ [GeV^]. Ordi- 
nate: St^-k+ [GeV^]. 
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Figure 4.7: Lines where = ±1, ±0.8, ±0.6, ±0.4, ±0.2, ±0.1. Abscissa: s^-^+ [GeV^] 

Ordinate: St,^k+ [GeV^]. 
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Chapter 5 

Guidelines for Dalitz plot analysis 



5.1 Surrogate scattering laboratories 

Experimentally, collisions with an unstable target (and beam) particle like irn vrvr and 
Ktt — ^ Kn are not feasible. Chew and Low fHI proposed that one can study such reactions 
in scattering of the beam particle off a stable target whereby the beam particle actually 
scatters with a virtual particle that is exchanged between the beam and the target particle. 
Two typical reactions of this type which serve as surrogate laboratory for scattering of 
unstable particles are n^p ir^iT^n and K"p — > K~TT'^n, see also section IB. 2. HI Another 
possibility in a similar vein is to study e. g. the scattering of two pions in the final state 
of a three-body decay such a.s J/ip ^ (pTm where as a good approximation the can be 
assumed not to interact with the pions [TTJ p. 1188, p. 1193]. As a by-product of the 
extensive CP-violation studies in B decays like K^n^n^ there is now a lot of new 

data available from a similar type of surrogate laboratory for final-state scattering of pions 
and kaons. One great advantage of using this new data is the huge amount of events such 
that it should be possible to interpret structures of the corresponding Dalitz plots that in 
earlier Dalitz plots would not be distinguishable from statistical fluctuations or not visible 
at all. 

Another feature of the reaction B^ K^n^n' is that contrary to the case of J/tp —>■ 
(pn-K the approximation that only two of the three final state particles interact is not 
applicable. Rather the pairwise interactions between all three final state particles are 
of the same order in strength. However, as an approximation one can assume that the 
interactions can be regarded as a superposition of three pairwise interactions and neglect 
the simultaneous interaction of all three particles in the final state, see equation ()5.H1 . 

1 will propose some guidelines for a Dalitz plot analysis for the three-body decay B^ 
K^TX^Ti~ . I concentrate on this example because it allows to show that an analysis in 
terms of a sum of Breit-Wigner terms and a non-resonant contribution a la Belle and 
Babar El naay be at odds with what is known about the involved two-body systems 
and insufficient to extract new or confirm old information about those. 

The goal of investigating the strong interactions of the final state particles of the three- 
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body decay defines a region of interest in the Dalitz plot. Controversial issues in the 
hadron spectroscopy of pions and kaons include the a, the k, the /o(980), the /o(1370), 
the /o(1500) and the i^*(1430). These are structures in the tttt and A'vr mass-region from 
zero to about 1500 MeV. These mass intervals corresponds to the lower left corner of the 
Dalitz plot of — * K~^n~^7r~ in the usual representations of e.g. Belle. ^ 

5.2 Three-body decay and two-body amplitudes 

As an approximation one can think of a three-body decay like -k^-k^K^ as the 

superposition of three amplitudes for the production of a two-body state with an accom- 
panying spectator, i. e. a particle which apart from its mere presence does not make any 
difference for the interactions taking place, 

+ T{B+ ^7i"{7r+K+)). (5.1) 

Each of the production amplitudes on the right-hand side can be related to the scattering 
amplitude of the respective two-body system. The relation that I will give is distinguished 
by that it is by construction consistent with unitarity constraints as discussed in the fol- 
lowing, cf. [niiHni. 

Here and in the following I use with refs. ED] T elements for discussing unitarity 
constraints. In my conventions the relevant quantity would rather he Ai. M. and T do, 
however, satisfy the same form of unitarity constraint, see equations (j3.96p and ()3.124|) . 
Therefore, the difference between T and M. is irrelevant here. 

5.2.1 Unitarity and production amplitudes 

To see how the production amplitudes of the right-hand side of equation ()5.1|) relates to 
the respective two-body scattering amplitudes, let us take as an example the amplitude 
T{B~^ —>■ IT' [ir'^ K~^)) . The unitarity constraint takes into consideration the channels B~^, 
71^ {n^ K^) and inelastic channels of the scattering of 7r~^K~^ (with 7t~ as spectator) like 
7r^(7r'*^/i''^7r^7r+). I enumerate the channels as 

l = fi+, 2 = 7r-(7r+i^'+), 3 = 7r-(7r+A'+7r-7r+), (5.2) 

where the dots stand for the prescription to enumerate all further inelastic 7r"'"i^'"'"-scattering 
channels. (There are no other channels with only two-particles to which 7t~^K~^ couples; 
no other two-particle state has the same quantum numbers as n~^K~^.) 
I use the unitarity constraint T'^T = i{T'^ — T) in the form 

\mT,, = ]^Y.^;,T,,. (5.3) 

k 

""^Note that Belle has s^-j^Jr as x-axis and s^-^+ as y-axis while in my plots I have (unintendcdly) 
reversed the assignment. 
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This means that I assume that the states can be properly discretized (cf. section 13.2.111 
and that time reversal invariance is given, ^ which I take to imply that the T matrix is 
symmetric, in other words 

= T*. (5.4) 

The meson consists of an up- and an anti-bottom quark. It has thus bottomness 1. 
Strong interactions conserve the flavors. Therefore, since the 5-meson is the lightest state 
with bottomness 1 it can only decay by electro-weak interactions. The elements Ti^ {k ^ 1) 
representing a transition from channel 1 (the B~^) to some other channel are therefore of 
order of magnitude of the electro-weak coupling constants, which in the energy region to 
be discussed here are small compared to the coupling constant of the strong interaction. 
In the following formulation of unitarity constraints I will neglect all terms quadratic in 
electro-weak couplings. To keep track of the different orders of magnitude of the electro- 
weak couplings I use a superscript 'w' for electro-weak amplitudes and a superscript 's' for 
strong amplitudes. 

The 5* element for interactions in channel 1 and transitions from channel 1 to any other 
channel is purely weak, 

Su = Sr,. te {1,2,3,...}. (5.5) 

On the other hand, interactions in and transitions between channels 2, 3, . . . are described 
by a Hamiltonian that is a sum of an electro-weak and a strong Hamiltonian, 

m^ = Ht^ + HZ, l,me{2,3,...}. (5.6) 

The corresponding S operator is therefore to leading order in Hf^ given by (see [221 P- 108]). 

S = lim exp{-iHim.{2t)) ^ lim (1 - ziy/;;t) exp{-iH^^{2t)). (5.7) 

t—*oo >oo 

For the T/^ {I, m G {2, 3, . . .}) elements we obtain 

Tim = -^iSim - Sim) ~ lim -tliSim - lHrm{2m5im + iT^m) ' ^Im] 

= Tlm^Tm + (5-8) 

with 

= lim(-i7-(2t)) and = Hm exp(-zi7- (2t)). (5.9) 

t—>00 t—fOO 

The factor gives a weak phase to the strong amplitude. This phase can be neglected in 
the present context as an overall phase. It can be absorbed by a redefinition of the state 
vectors. Not keeping track of this overall phase and neglecting quadratic and higher order 
terms of the weak T elements we then have in and between channels 2, 3, ... 

Tji = T;"^+T^'„ j,te{2,3,...}. (5.10) 

^Because of CPT invariance the assumption of time reversal invariance is incompatible with the possi- 
bility of CP violation. For the issues discussed here the assumption of T invariance is simplifying but not 
necessary. CP violating phases can be neglected in the present context as overall phases. They can be 
absorbed by a redefinition of the state vectors. 
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The unitarity condition for the production amphtudes, T^*^ [k ^ 1), reads 

2imTr, = J]ri;r,, (5.11) 

r=l 

= TiTT- + J2 Tir* (T4 + r,) (5.12) 

r=2 

= X^TrT4, (5.13) 

r=2 

where in the last step terms quadratic in the weak couphng have again been dropped. 
If the production amphtudes are written as, 

T,, = Tr, = J2^nT:,, Ml, (5.14) 

n=2 

where, importantly, the a's are real, the constraint of equation (j5.13|l is satisfied by con- 
struction; provided that the strong amplitudes T^^ satisfy the unitartity condition among 
themselves, i.e. 

2ImT;^ = ^T-T,^^, /c,/,me{2,3,...}. (5.15) 

k 

Indeed we then have: 



2ImT- = 2Im ^ a„7^4 = «n2Im7^4 = E E ^-^'^^^ 

n=2 n=2 n=2 r=2 

r=2 n=2 r=2 



5.2.2 Elastic region for tt tt'^ and tt 

Because of equation ()4.16|) not all three two-body cms energies can be low. The region of 
the Dahtz plot we are interested in is e. g. characterized by low values of S7r-7r+ and s-^-k+ 
and high values of St,+k+, roughly 

GeV< v/i~;T< 1.6 GeV, (5.17) 
GeV < <1.6 GeV, (5.18) 

3 GeV < v^i^T^ < 5 GeV. (5.19) 

From analyses of vrvr scattering (e.g. CERN-Munich (211, see figure lower panel) it 
is known that indeed the inelasticity is almost 1 up to 1.5 GeV, except for a small energy 
region around 980 MeV. This dip in the inelasticity may be due to the decay of /o(980) 
into TTTT as well as into KK. 1 neglect this inelasticity in the hope that the resulting ansatz 
is while simpler still able to represent the most relevant features of the process. 
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For the production amplitudes T{B^ (tt^tt and T{B^ {n^ K^)7r'^) we can 
therefore use an elastic unitarity condition that takes only elastic strong amplitudes to be 
non-zero. The respective production amplitudes have then the simple form Ti^ = T^*^ = 
afcT4, i.e. 

T{B+ ^ (7r-7r+)Ji+) = a^-^+T%n~n+ ^ 7r-7r+) 

T{B+ (7r~A'+)7r+) = a^- k+T' {tx~ 7r~A'+). ^ ' ^ 

This form of the production amplitude shows that, since ak is real, the production am- 
plitude and the strong amplitude of the interaction of the final state particles have the 
same phase. This result is known as Watson's theorem, final state interaction theorem 
or just elastic unitarity. The name "final state interaction theorem" is justified by that 
the particular form of the production amplitude (eq. ()5.14j) ). which is a solution for the 
unitarity constraint can be interpreted to the effect that the decay from channel 1 into a 
channel /c 7^ 1 is a sum of amplitudes representing the decay from channel 1 into channel 
n ^ 1 with coupling strength a followed by a strong final state interaction represented by 
l3Uj. In the original Watson theorem J3] it is explicitly two potentials that are 
considered: The potential responsible for the production of the hadrons, and the potential 
of the interaction of these hadrons in the final state. 



Singular couplings 

Since the partial wave amplitudes satisfy a unitarity condition of the same form as the 
elements of the T matrix, see equation ()3.109|) unitarity for production amplitudes re- 
quires also each partial wave and isospin component to be of the form of equation ()5.14|1 . 
Thus we obtain unitarity constraints for the elastic amplitude in terms of phase shifts (cf. 
equation ()3.114j) ). for example: 



M-{{B+ {tx-tx+)K+) = a-i{Mi 



■^^"'tt 7r+ ^ tt tt^ 



(5.21) 



As discussed and emphasized in [HI p. 1188] and [HH] the a's are not necessarily regular 
functions. They may be such that zeros of the elastic amplitudes are removed. This is the 
case if a behaves like 1/ sin 5 in the region where 5 is a multiple of vr. Also the a's may 
introduce new zeros. It is adequate to redefine the real coupling constants 

2a/ sin 5/ = a/, (5.22) 

such that 

X/(S+ ^ (7r-7r+)i^+) = a^jc'^'. (5.23) 

In ref. |1H p. 1189] it is emphasized that, while the production and the scattering 
amplitude do not necessarily have the same zeros, they do indeed have the same resonance 
poles. Resonance poles are in this sense universal: "A further consequence of unitarity. 
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vital for our discussion, is that it requires that resonance poles be universal, i.e. , a given 
resonance pole occurs at the same complex energy Eji in all processes to which it couples. 
This is automatically built in our solution, [...]" 

I share the view that the universality of resonance poles is indeed an important feature 
of the solution given in ref . jTI] . The universality gives the resonances an identity. However, 
since poles are complex quantities they are not observable and their determination therefore 
requires an extrapolation from real to complex values of energy. This in turn presupposes 
a lot of theoretical input. 

More specifically, I do not see how the universality of the resonance poles is implied by 
unitarity. With the solution of ref. [TT] the scattering and the production amplitude do 
have the same poles; it respects the universality of the resonance poles. But I cannot see 
an argument why this is the only way to satisfy the unitarity constraint. 

5.2.3 High energy amplitude for i^+7r+ 

In the region of interest defined in equation ()5.17|1 we have to consider the high energy 
behavior for the amplitude T{B'^ tt~ {7t~^ K~^)) in contrast to the amplitudes T{B^ 
{n~'n'~^)K~^) and T{B~^ {tt~ K~^)7r'^), where we are concerned with low cms energies. 
Unitarity as formulated in section l5.2.1l is satisfied by the following ansatz (see eq. ()5.14|) ). 

+ a^+K+n-n+T'{7r+K+'K-n+ 7r+K+) + ... (5.24) 

No resonances 

The two particle state tt'^K'^ has isospin / = 3/2 and has electric charge 2. A qq resonance 
cannot have this electric charge and this isospin. Baryon resonances can, e.g. A"*""*", but 
because of baryon number conservation 7r~^K~^ baryon is not possible. To conserve 
baryon number, baryonic resonances should be produced as pairs of baryon and antibaryon. 
Also there is the a priori possibility of resonance of four quarks. In any case, no resonances 
in the 7r~^K~^ channel are known as of today. 

The S-wave phase shift with I = 3/2 is at least in the interval 0.8 GeV to 2 GeV 
negative [21]. This supports the claim that there are no resonances in this channel; a 
repulsive force cannot lead to resonances. 

Forward peak in Ktt scattering? 

I do not see any reason why the ir'^K^ total cross section should vanish. By the optical 
theorem (see section 13. 9p then the imaginary part of the elastic amplitude in the forward 
direction is non-zero and with it the absolute value of the elastic amplitude. So although 
there are no resonances in the 7r~^K~^ channel we have reasons to expect that the elastic 
amplitude is not zero in the region of the Dalitz plot under consideration. Contrary to 
this expectation the amplitude in the 7r~^K~^ channel is set to zero in the default-model of 
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ref. jHj, see section IH.H. 21 There is in principle the possibility that a7r-(7r+/<+) tends to zero 
in such a manner that the product Ut^- i^t^+ k+)T [ti^ -k^K^) is zero at values for the 
n^K^ invariant mass of about the B mass. However, since I know of no reason why this 
should happen to be so I assume that also the product is different from zero. 

Not only we have reasons to expect that the K^tt^ amplitude is not zero. There is also 
experimental indication that the K^tt^ amplitude as defined in chapterElis forward peaked 
at energies relevant for the Dalitz plot analysis of B^ 7r~TT^K~^: The Lass collaboration 
obtains in ref. J7] results about K^ir'^ scattering by studying the reaction 

K~p A'^TT+n. (5.25) 

Basically, the method to extract Ktt scattering data from such reactions is the one of Chew 
and Low One selects the events with low momentum transfer t. In these events the 
exchange of a virtual pion in the t channel is the dominant contribution to the interaction 
of the proton and the kaon. By an extrapolation from virtual t's to the pion pole one can 
then obtain information about the scattering of pions and kaons. The scattering angle 6 for 
Kir — > Kn is essentially the Gottfried- Jackson angle ^gj of ths reaction K~p — > K~Ti^n. 

Figure 15.11 shows a distinct enhancement of the K'n^ cross sections in the forward 
direction (cos^gj = 1)^ and energies larger than 2 GeV. This is particularly so for events 
with high t (bottom panel) where the exchange of heavier particles as the p dominate over 
the pion contribution. But also in the case of relevance here — low t's, pion exchange — the 
events at energies around 2.5 GeV and higher are concentrated at values of 0.5 < cos6'gj < 
1. 

The A'+TT"*" scattering is purely isospin / = 3/2. The K^ti^ amplitude, studied in [T7] . 
also contains an J = 3/2 component. It is not clear to me whether the forward peak comes 
only from the isospin components with J 7^ 3/2. As a hypothesis I nevertheless assume 
the amplitude T^{ti^K^ ti^K^) at energies of about 4-5 GeV, which are of interest 
here, to share the feature of being forward peaked (or backwards, see footnote El) . From 
these rough arguments I dare conclude that as a hypothesis we may write for the tt^K^ 
production amplitude a non-resonant amplitude with an enhancement of its absolute value 
in the region where |2;7r+K+| ~ 1. This is near the boundary of the Dalitz plot. In particular 
for pair masses of about 4-5 GeV for ir^K^ the enhancement would be near the lower left 
corner of the Dalitz plot. I denote the amplitude by 

T(5+ ^ 7r-(7r+A^+)) = (5.26) 
5.2.4 Pronounced S-wave? 

Compared to the K and the D the B has a much higher mass (m/^ ^ 495MeV, m£, ^ 
1900MeV, rriB ~ 5300MeV). The phase space of the three-body final state is there- 
'll am not sure if 9gj corresponds indeed to the forward scattering angle or rather to the one backwards. 
However, what is of importance here is only whether the K~^Tr~^ amplitude can generate an enhancement 
near the boundary of the Dalitz plot for B+ 7r~7r+isr+; and both, 9 = and 9 = 180 degrees are lines 
on the boundary, see figures 14.81 
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Figure 5.1: The cosine of the Gottfried- Jackson angle of p ^ K n^n versus the pair 
mass oiK-Ti+. Top paneh \t\ < 0.2 GeV^ Bottom paneh 0.2 < \t\ < 1.0 GeV^; from [Ilj. 
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fore much larger in three-body decays of the latter. The region of interest (see equa- 
tions ()5.17j) ff.) for present problems of hadron spectroscopy is given by relatively low 
masses of the 7r~7r+ and tt~ pair. By equation ()4.18|) the pair mass of tt^K^ is there- 
fore about 5 GeV. If the tt^K^ amplitude is indeed forward (or backward) peaked as 
suggested in section I5.2.3| this could lead to an enhancement near the boundary of the 
Dalitz plot in the lower left corner. I would be interested in seeing whether such an en- 
hancement explains at least in part the purported observation of an important S wave 
signal. The risk of mistaking the tt^K^ forward peak for a strong scalar signal of e. g. 
/o(980) is given because the bands or dips expected for such a state is very close to the 
left boundary, see figure l5?2l 

5.2.5 Resulting ansatz 

Putting together equations ()5.H) . ()5.2m) and ()5.26p we obtain for the three-body decay 
amplitude 



TT TT 



+ a^-K^T^{Tx-K^ ^ 7r-A'+) + A^'J^t'^ (5.27) 

Performing isospin decomposition (with conventions of the Particle Data Group) we get 

T{B+ 7r-TT+K+) = /6T2(7r7r) + l/2Ti(7r7r) + l/3To(7r7r)] 

+ a^-KA^/3Ts/2{7c-K+) + 2/3Ti/2(7r-A'+)] (5.28) 

^ ^boundary ^^+^^+^_ 

Partial-wave decomposition for the tt^tt^ and the 7r~K~^ amplitudes (taking into ac- 
count only S and P wave, i. e. J = and J = 1) further yields^ 

T(5+ ^ 7r-7r+A'+) = -^(a^-^+[l/6r°(7r7r) + 1/2 cos^^-^+ri^(7r7r) + l/3r°(7r7r)] 

+ a^-,,+ [l/3{T3%(7r"A'+) + cos^^-,,+T3%(7r"A'+)} 29) 

^ ^boundary ^^+^+^_ 

Because the (real) coupling constants a may be singular (see section I5.2.2|) it is appro- 
priate to define the following (still real) coupling constants: 

(a^-^+yj = 2a^-^+ sin(5/(7r"7r+), 

7 7 _ -i_ (5.30) 

(Q!^-i^+)j = 2a^-i^+ sin(5j (tt K~^). 



"^Because of Bose symmetry there is no component with J = and / = 1 or J = 1 and / = in the tttt 
channel. 
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Then equation ()5.29|1 reads 

T{B+ n-n+K+) = -L (i/6(a,_^+)0e''^°(--) + 1/2 cos e,-,+ (a,-,+ )Je*'^^("") 

47r 

+ l/3{(a..^+)°/2e^'3/.(--^+) + cose..x+(a.-^+)^/2e^'3/.(--^+)} 
+ 2/3{(a.-^+)?/2e</^(^"^^) + cos^.-^+(a.-A^+)}/2e^'^V^(^^^^)}) 

(5.31) 



In this ansatz the phases of the n~7r~^ and tt~K~^ amphtudes are not free fit parameters 
but are the strong phases of the respective scattering amphtude. The phase Sq^hit) for 
example is shown in figure (top panel). In the Dalitz plot event distribution, which is 
essentially given by |Tp the phases show up in interference terms of type 



cos(5^(7r7r) -(5°(7r7r)). (5.32) 



5.3 Two examples 



In the following two subsections I briefly report two analyses of K^-k'^ti' and 

— >■ K~^K~^K~ . The idea is to give an impression of how different approaches lead to 
different conclusions. 

5.3.1 Minkowski and Ochs 

In ref. jHl table VII] the branching fraction for the decay B^ — > K'^7r~^7r~ via the resonance 
/o(980) is found to be 7.55 x 10"^ while in Table 6] it is concluded to be 19.3 x 10"^ 
One of the reasons for such discrepancies in the analysis of the three body decays B^ 
K^7T~^7T~ and B'^ K^K^K~ is the role played by the object called "non-resonant 
background" or "gb" respectively. In [0] , the gb interferes in the ti^tt' channel destructively 
with both /o(980) and /o(1500) such that the peaks in the corresponding mass projection 
of the Dalitz plot appears smaller than might be expected or even as a dip, see [0] for 
details, in particular figure 4. In the channel K~^K~ , on the other hand, the gb interferes 
constructively with the /o(1500) such that, contrary to the expectations, the number of 
events in the region of the /o(1500) in the K^K~ mass projection is actually much higher 
than in the tt+tt" projection, see [Oj, in particular figure 5. 



Large background phase, no /o(1370) 

Minkowski and Ochs fit only the mass projections on ti'^ti" and K^K~ . According to 
them a fit that takes into account the two-dimensional Dalitz plot distribution would be 
preferable but as a first exploration a fit of the one-dimensional mass projection should give 
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a rough picture of the dominant aspects. The ansatz for the rate per pair mass interval is 
roughly the following (m = y/s), see jHl eq. 12] 



dm 



hh 



hh = 7T+7r-,K+K-. (5.33) 



bg 



This is a superposition of an amplitude representing the broad resonant glueball, |Tgb|e 
and two amplitudes representing the resonances /o(980) and /o(1500). The two resonances 
are represented by Breit-Wigner terms^ with an energy dependent width which, moreover, 
is characterized by a particular shape function G, 

Tfi r 

Ta = — . Vn ( \V « = ^?^/o(980),/o(1500). (5.34) 

- - imaTa[l + Ga("^)) 

The relative phases between the three resonant contribution are not free fit parame- 
ters but share the common factor of e^^''^ . This has as a consequence that the /o(980) and 
/o(1500) are not typical Breit-Wigner resonances the phases of which move rapidly through 
7r/2 but are rotated in the Argand diagram by two times the phase of the glueball ampli- 
tude at resonance (7r/2). This ansatz is motivated by the interpretation of the isospin S 
wave as the red dragon, the broad 0"^^ glueball interfering destructively with the narrower 
resonances /o(980) and /o(1500), see figures 1^75^13. 71 

Interestingly enough, the superposition of all three amplitudes leads to an enhancement 
in the region of m = ^Jshh = 1300MeV, which the Belle Collaboration (see section I5.3.2|l 
would rather interpret clS db SI gnal from the controversial /o(1370). 



5.3.2 Belle Collaboration 

I do not discuss here the event reconstruction, background suppression, and corrections 
for efficiency and resolution of ref. [S|.^ I only report briefly the ansatz for fltting of the 
signal. The model that is used to determine the branching fraction in question consists of 
the following ansatz for the S-matrix element: 

Saj{K^7i^7^~) = aK*e'^''*Aii7T+K+'K-\K*{892)^) 

+ e'^"o A (vr+K+TT- 1 K* ( 1430) ) 
+ ape'^'^Ai{K+TT+TT~\p{770f) 

+ af,e'^foAo{K+7i+7i-\fo{980)) (5.35) 

+ af,e''fxAj{K+7i+n~\fx) 

+ a^^,e'^^r^oAo{K+7i+7i-\xco) 
+ Ar(/^+vr+7r-), 

^The T of ref. jH| is normalized differently than my A4. Cf. e.g. the expression for a Brcit-Wigncr 
resonance, equation H3. 118(1 . 

^An important point to be checked would be how the pairs of signal events are formed for different 
mass projections, as discussed in section I?. 1.21 
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with [HI eq. (9)] 

Anr{K+n+n-) = ^■•e-^^^e*^!" + a^'e-'''^'e'^^^\ (5.36) 

Here I adopt the notation of ref. |H]. In this notation the order of the symbols ^K'^\ 
'vr^' and 'vr"' in the brackets of A indicates that the first two terms in equation ()5.35|) 
denote an amphtude in the K~^7r~ channel whereas the second to sixth term denote each 
an amplitude in the 7r~^7r~ channel. As to the non-resonant amplitude Am, the subscripts 
'13' and '23' show that this amplitude is the (coherent) sum of one amplitude in the K^7r~ 
channel and one amplitude in the tt+tt" channel. 

In order to estimate model dependent uncertainties a fit was also carried out using a 
non-resonant amplitude with a contribution 

^nr^-as^,^^Sr (5.37) 

from the K~^7r~^ channel. In the default model, called Knir — Co, to determine the branching 
fractions, however such a contribution is not present. 

Let me spell out more explicitly what the ^'s are. I take as an example ^i(7r"'~/r^7r~|i^'*(892)'^). 



A(7r+A'+7r"|A'*(892)°) oc 



A/2 _ - lMRr]^{si3) 

{Ml -ml) {ml -ml] 



X Si2 - S23 + 



■§13 



(5.38) 



with 



r^(si3) oc — —= ^ ^ ^ . (5.39) 



5l3 



The first factor on the right-hand side of equation ()5.38p is a relativistic Breit-Wigner 
function with energy-dependent width r}j(si3), the second factor describes the angular 
distribution of the three decay products for a resonance with total spin J = 1. 

For a J = resonance, like for example the /o(980), the amplitude A takes a simpler 
form. The factor for the angular distribution is 1. For Mr, ^/s23 ^ m^^ the energy depen- 
dence of the width can be neglected, r5j.(s23) ~ Tr. The corresponding density of events 
in the Dalitz plot is proportional to A"^, which gives in case of no interference a band of 
width r|. parallel to sis the and perpendicular to the S23 axis. 



Non-resonant background and /x(1300) 

In contrast to the ansatz of Minkowski and Ochs the analysis by the Belle collaboration 
lacks a clear interpretation for the introduction and parametrization of their non-resonant 
amplitude A^r- As to the controversial /o(1370), the Belle collaboration finds in their 
default model Knir — Co that the "the mass and width of the /x(1300) state obtained 
from the fit with the model Kirir — Co are consistent with those for the /o(1370)" |H, p. 15]. 
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5.4 Conclusions 

In conclusion, I propose the following guidelines to begin with a Dalitz plot analysis of the 
decay ti'ti^K^ or comparable. It goes without saying that I do not mean to have 

derived strict rules that are necessary and sufficient for a successful analysis. I only try to 
stress some points that might be neglected in common fit procedures and want to suggest 
some steps that seem appropriate in the light of the issues discussed in precedent sections 
and chapters. 

5.4.1 Assumptions to be checked 

The guidelines I propose rest in particular on the following assumptions, which remain to 
be checked. 

• The two-body angles defined in section the Gottfried- Jacks on angle and the angle 
9 in two-body cms scattering are directly related. 

• For the tt^tt^ and ir^K^ channel one can use as a constraint elastic unitarity. In- 
elastic channels can be neglected. 

• The 7r+A'+ amplitude has a forward (or backward or both) peak. 

• The coupling of of the B to the tt'^K'^ does not vanish, see page ESI 

5.4.2 Guidelines 

Forming triplets of signal events. One must not only held fix the signal to background 
ratio but form triplets of signal events for the entire analysis including in particular different 
mass projections, as discussed in section 1^.1.21 

Two-dimensional binning. Dalitz plot analysis of a decay like -k^-k^K^ differ 

from older Dalitz plot analysis in that they have not the principal aim to determine spin 
and parity of the decaying particle from characteristic patterns in the plot, see figures fTTTl 
and 11.21 Rather the aim is to fit the "landscape" over the plot. To this end two dimensional 
binning of the data and representation as a lego plot is more appropriate than a distribution 
of points. 

Do not fit whole Dalitz plot. The interior of the Dalitz plot boundary is approximately 
3 times larger than the actual region of interest in the lower left corner with St^-t^+ and St^-k+ 
between zero and 1.6 GeV, roughly. A fit of the hadronic structures under consideration 
should only be done in this region of interest. If one fits the whole Dalitz-plot, a certain 
set of parameters may score better in a test, but not because it represents better the 
structures in question but the huge rest of the plot, which is irrelevant for the present 
problems of hadron spectroscopy. 
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Elastic phases consistent with scattering data. Because of the at least approximate 
validity of elastic unitarity constraints (Watson's theorem, see section I5.2.2|l the phases of 
the amplitudes T{B~^ {7i^7i'^)K~^) and T{B^ {t\:~ K^)ti^) should be (approximately) 
equal to the phase shifts established in tttt and Kn scattering analyses. 

Forward (or backward) peak from K~^tz^1 Without good reasons the amplitude 
T{B^ -K^^-K^K^)) should not be set to zero. Quite the contrary may be true: This 
amplitude may be non-vanishing and generate an appreciable enhancement at the lower 
left boundary of the Dalitz plot. If such a signal is indeed present, care should be taken to 
separate it from searched for S wave signals; see section 15.2.41 

K*{S92)^ and Zt^-j^+ = 0. The uncontroversial ii'*(892)° resonance can serve as refer- 
ence for fixing overall coupling strength and phases. As a spin 1 particle the contribution 
from this resonance should vanish with 2;^-/^+. 

Peak from i^+Tr"'" and zero from 1^*(892)°. I could not establish in detail where the 
line with ZTr-K+ passes through the lower left corner of the Dalitz plot, see figure I^Tl Since 
the signal from i^*(892)° should vanish there the non- vanishing of the K^tt^ amplitude 
could most be noted. 

Charge and flavor conjugated reactions. Analyses of charge and flavor conjugated 
reactions to -k^-k^K^ like 

S+, 5", 5°, fio ^ tttt/s:, r]r]K, 7]ttK (5.40) 

may give additional input to decide between alternative fits, and contributions from dif- 
ferent isospin components. 
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